ASYMPTOTIC GLUING OF ASYMPTOTICALLY HYPERBOLIC SOLUTIONS 
TO THE EINSTEIN CONSTRAINT EQUATIONS 
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Abstract. We show that asymptotically hyperbolic solutions of the Einstein constraint equations 
with constant mean curvature can be glued in such a way that their asymptotic regions are con- 
nected. 



1. Introduction 

One of the most useful ways to produce new solutions of the Einstein constraint equations is via 
gluing techniques. The standard gluing construction is the following: We presume that (M, g, K) 
is an Einstein initial data set, with M a smooth n-dimensional manifold, g a Riemannian metric 
on M, and K a symmetric tensor field on M. We further assume that this set of data satisfies the 
(vacuum) Einstein constraint equations 

(1) dw g K - V Tr g K = 0, 

(2) R(g) - \K\ 2 g + (TvgK) 2 = 0, 

which are the necessary and sufficient conditions for (M, g, K) to generate a spacetime solution of 
the (vacuum) Einstein gravitational field equations via the Cauchy problem [6]. (Here div 9 is the 
divergence operator, Tr g is the trace operator, R(g) is the scalar curvature, and | • \ g is the tensor 
norm, all corresponding to the metric g.) 

Choosing a pair of points pi,P2 £ M, one shows that there is a family of new solutions 
(M e , g e , K E ) of the constraint equations in which (i) M e is obtained from M by (connected sum) 
surgery joining p\ and p2 and (ii) outside of a neighborhood of the connected sum bridge in M e , 
the data (g £ , K £ ) can be made as close as desired to (g, K) (in a sense to be made precise later) by 
taking e sufficiently small. We note that this gluing construction allows for the possibility that the 
manifold M consists of two disconnected components; then if p\ is chosen to lie in one of the com- 
ponents and p2 in the other, the new glued solution effectively connects two disconnected solutions 
of the constraint equations. 

The mathematics and the utility of the gluing of solutions of the Einstein constraint equations 
are discussed in a series of papers [10, 14-16], which show that gluing can be carried out for 
a wide variety of initial data sets: they can be compact, asymptotically Euclidean ("AE"), or 
asymptotically hyperbolic ("AH"), and they can be vacuum solutions or non- vacuum solutions 
with various coupled matter fields. This past work shows that in some cases the gluing can be 
done so that the glued solution exactly matches the original one outside the gluing region, so long 
as certain nondegeneracy conditions ("no KIDS") hold at the points of gluing. When this can be 
done, the gluing is said to be localized. More generally, the glued solutions may not exactly match 
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the original ones outside the gluing region, but can be constructed so that the data set is arbitrarily 
close to the original solution away from this region; this is called non-localized gluing. 

Say one chooses a pair of (disjoint) asymptotically hyperbolic solutions of the constraints and 
glues them at a pair of points satisfying the necessary conditions, as described in either [15] or [10]. 
If each of the original AH data sets has a single (connected) asymptotic region (as described below 
in Section T2.ip . then the glued data set, which is also asymptotically hyperbolic, necessarily has two 
disjoint asymptotic regions. If we are working with AH initial data sets which are viewed as data 
on partial Cauchy surfaces that intersect null infinity in an asymptotically simple spacetime [22], 
then the existence of multiple asymptotic regions is problematic for physical modeling. 

In the present paper, we show that one can glue asymptotically hyperbolic solutions of the 
constraint equations in such a way that in fact the asymptotic region of the glued data set is 
connected. The idea, which is modeled after the studies of Mazzeo and Pacard on gluing asymp- 
totically hyperbolic Einstein manifolds [21], is to use the conformally compactified representation 
of asymptotically hyperbolic geometries, which models a complete, asymptotically hyperbolic man- 
ifold as the interior of a compact manifold with boundary (see Section [2.11 below). The boundary 
of this manifold, which we call the ideal boundary, is not part of the physical initial manifold, but 
represents asymptotic directions at infinity. The gluing is done using points p\ and p2 lying on 
the ideal boundary. In this context, we will use the term asymptotic region to refer to any open 
collar neighborhood of the boundary, with the boundary itself deleted. If the original manifold has 
a single connected asymptotic region, then so does the glued manifold. (See Fig. [U) If the original 




Figure 1. Figure 2. 

manifold has two disjoint connected asymptotic regions with one boundary point chosen on the 
ideal boundary of each region, then the glued manifold will have a single connected asymptotic 
region (Fig. [2]). (Although the ideal boundary in Fig. Q] appears to be disconnected, in the cases 
of interest the boundary of the glued manifold will be a connected sum of connected 2-manifolds, 
which is always connected.) 

The results we present here do not hold for general AH initial data sets. We require that the 
data have constant mean curvature ("CMC"), in the sense that Tr 5 K is constant on M. Also, 
our results thus far provide sufficient conditions for (asymptotic) non-localized gluing. In a future 
paper, we hope to both eliminate the CMC restriction, and find conditions which are sufficient for 
the gluing to be localized. 

To set up our work here, we start in Section [2] with a definition and discussion of asymptotically 
hyperbolic initial data and their polyhomogeneous behavior in the asymptotic region. The section 
continues with a brief description of the conformal method for generating solutions of the constraint 
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equations and the simplifications of the method which occur for constant mean curvature data. The 
conformal method is discussed here because it plays an important role both in constructing basic 
examples of AH initial data [1] and in carrying out our gluing procedure. We conclude Section 
[2] with an overview of those aspects of [1] on which our work relies. We state our main theorem 
in Section [3J We then carry out the first part of the gluing construction ( "splicing" ) in Section 
HI producing a family of initial data sets depending on a small parameter e, which satisfy the 
constraint equations approximately. Our splicing construction is modeled after the construction of 
asymptotically hyperbolic Einstein metrics in [21]. 

The proof of the main theorem relies on the use of weighted Holder spaces. We define and 
discuss these spaces in Section following [18]. In the rest of Section we study certain elliptic 
operators which act on weighted Holder spaces, and prove that they are invertible with norms of 
their inverses bounded uniformly in e. One of the key steps here is a blow-up analysis argument 
similar to that of [21]. Our analytical results are applied in Section [6] to correcting the traceless part 
of the glued second fundamental form. Finally, in Section [7] we use results of [12] and a contraction 
mapping argument to solve the Lichnerowicz equation and complete the proof of the main theorem. 
Section [8] contains concluding remarks. 

2. Preliminaries 

In this section we define and discuss examples of asymptotically hyperbolic initial data sets. We 
also review the conformal method for creating solutions of the constraint equations. 

2.1. Asymptotically Hyperbolic Initial Data Sets. The model for asymptotically hyperbolic 
initial data sets is the hyperboloid in (3 + l)-dimensional Minkowski space, 

M = {(x°, x\ x 2 , x 3 ) : x° > and (x ) 2 = (x 1 ) 2 + (x 2 ) 2 + (x 3 ) 2 + 1}, 

together with its induced Riemannian metric g (which is a model for the hyperbolic metric of 
constant sectional curvature —1) and its extrinsic curvature K = g; it is straightforward to check 
that (M, g, K) satisfies {TJ and ©. 

Roughly speaking, an asymptotically hyperbolic initial data set (M, g, K) is one in which the 
Riemannian metric g is complete and approaches constant negative curvature as one approaches 
the ends of the manifold, and the extrinsic curvature K asymptotically approaches a pure trace 
tensor field that is a constant multiple of the metric. 

A more precise definition of asymptotic hyperbolicity is motivated by the Poincare disk model 
of hyperbolic space. In that model, hyperbolic space is given by the (smooth) metric 

(3) g= (l-|x| 2 ) 2 {{dxXf + ' ' ' + ' 

on the open unit ball. If we define the function p = i(l — |x| 2 ), then g = p 2 g is the Euclidean metric 
and thus extends smoothly to the closure of the open ball. Note that the function p is smooth on 
the closed ball, it picks out the boundary of the ball since p _1 (0) is equal to this boundary, and it 
satisfies the derivative condition \dp\g = \^d(l — \x\ 2 ) | = 1 on the boundary of the ball. 

With this example in mind, we make the following definitions. We suppose throughout this 
paper that M is the interior of a smooth, compact manifold with boundary M. A defining function 
for M is a nonnegative real- valued function p : M — ► R of class at least C 1 such that /9 -1 (0) = dM 
and dp does not vanish on dM. Given a nonnegative integer k and a real number a £ [0,1], a 
smooth Riemannian metric g on M is said to be conformally compact of class C k (or C k,a , or C°°) 
if there exists a smooth defining function p and a Riemannian metric g on M of class C k , C k,a , or 
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C°°, respectively, with g = p~ 2 g~\M- Smoothly conformally compact means the same as conformally 
compact of class C°°. 

If g is conformally compact and if in addition \dp\g = 1 on dM, then g is said to be asymptotically 
hyperbolic (of class C k , C k,a , or C°°, as appropriate). One verifies easily that any asymptotically 
hyperbolic metric of class at least C 2 has sectional curvatures approaching — 1 near dM (see [20]), 
and so indeed has the intuitive properties mentioned above. The boundary dM is called the ideal 
boundary, and dM together with its induced metric t*g (where t : dM ^ M is inclusion) , is called 
conformal infinity. Note that for a given AH metric, the choice of the defining function is not 
unique, so the geometry of the conformal infinity is only defined up to a conformal factor. 

Naively, one might hope to work with initial data sets (M, g, K) in which (M, g) is a smoothly 
asymptotically hyperbolic Riemannian manifold, and K is a symmetric 2-tensor field such that 
p 2 K has a smooth extension to M and is equal to a constant multiple of the metric on the ideal 
boundary. Unfortunately, however, it is shown in [2] that there are obstructions to finding solutions 
with this degree of smoothness, marked by the presence of log terms in the asymptotic expansions 
of g and K near the ideal boundary. For this reason, instead of smoothness we have to settle for a 
slightly weaker notion called polyhomogeneity, which we now define. 

A function /: M —* R is said to be poly homogeneous (cf. [20]) if it is smooth in M, and there 
exist a sequence of real numbers s, f +oo, a sequence of nonnegative integers and functions 
Jij e C°°(M) such that 

OO Qi 

(4) /~£X> ai 0ogp)% 

i=\ j=0 

in the sense that for any positive integer K, there exists a positive integer N such that the difference 

N qi 
i=l j=0 

is 0(p K ) as p — > 0, and remains 0{p K ) after being differentiated any number of times by smooth 
vector fields on M that are tangent to dM. It is easy to check that sums and products of polyhomo- 
geneous functions are polyhomogeneous, as are quotients of poly homogeneous functions provided 
that the denominator has no log terms with its lowest power of p (i.e., q\ = 0) and provided that its 
leading term / 10 does not vanish on the ideal boundary. A tensor field on M is said to be polyho- 
mogeneous if it is smooth on M and its component functions are polyhomogeneous in some smooth 
coordinate chart in a neighborhood of every ideal boundary point. We define a polyhomogeneous 
asymptotically hyperbolic Riemannian metric on M to be a polyhomogeneous Riemannian metric 
g which is also conformally compact of class at least C 2 . 
Now we come to the main definition of this section. 

Definition 1. A polyhomogeneous asymptotically hyperbolic initial data set (sometimes called a 
hyperboloidal initial data set) is a triple (M,g,K), in which 

(i) M is the interior of a smooth, compact manifold with boundary M ; 

(ii) g is a polyhomogeneous AH Riemannian metric on M; 

(iii) K is a polyhomogeneous symmetric covariant 2-tensor field on M with the property that for 
any smooth defining function p, p 2 K has a C 2 extension to M whose restriction to dM is 
a constant multiple of (the extension of) p 2 g there; 

(iv) the Einstein constraint equations (JT]) and fl2J) are satisfied. 
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This definition implies that g and K can be written in the form 

(5) g = p~ 2 g, 

(6) K = -g + p- 1 ]!, 

n 

where g is a poly homogeneous Riemannian metric on M that is of class at least C 2 (and thus has 
log terms, if any, only with powers of p greater than 2); r is a poly homogeneous C 2 scalar function 
on M whose restriction to dM is constant; and 71 is a polyhomogeneous symmetric 2-tensor field 
on M that is trace-free with respect to ~g and has log terms only with powers of p greater than 
1. A polyhomogeneous AH initial data set is defined to be CMC if the mean curvature function 
r = Tr g K is constant. We note that for such data r must be identically equal to 3. Indeed, the 
definition of CMC asymptotically hyperbolic data guarantees that the scalar curvature approaches 
the constant —n(n — 1) = —6 at the ideal boundary. Since the g-norm of p~ l ~p approaches zero at 
the ideal boundary, inserting ([6]) into ([2]) and evaluating in the limit at the ideal boundary implies 
that the mean curvature r must be identically equal to 3. Thus ([6]) becomes 

K = g + p = p~ 2 g + p _1 /Z. 

It is shown by Andersson and Chrusciel in [1] that CMC polyhomogeneous AH initial data sets 
exist in abundance. 



2.2. The Conformal Method for Finding Solutions of the Constraint Equations. Both 
the construction of [1] and our gluing construction here are based on the conformal method, which 
(along with the closely related conformal thin sandwich method) is the most widely used method 
for producing solutions of the constraint equation. We proceed by reviewing this method. 

We start by introducing two auxiliary differential operators which are involved in the conformal 
method. The first of the two is the conformal Killing operator V x which acts on vector fields X as 
follows: 

(7) (V x X) cd := l -C x X cd - i(div A X)A cd = \{V c X d + V d X c ) - ^V a X a X cd . 

The image of T>\ is contained in the space of symmetric 2-tensors which are traceless with respect 
to A. The formal adjoint of T>\ is 

(8) V* X T = -(dw x Tf; 

here (and throughout the paper) the symbol ' refers to raising an index. Another auxiliary operator 
we use is the elliptic, formally self-adjoint vector Laplacian operator 

L\X = (V* x o V X )X = - (div A [V x X)f . 

The conformal method for 3-dimensional manifolds (we restrict to n = 3 for convenience) is 
based on the Lichnerowicz- York decomposition of data [7] 

(9) g a b = 4> 4 Kb, 

(10) K cd = V>- 2 (y cd + 2(V x W) cd ) + i V 4 r X cd , 

where \ ao is a Riemannian metric, v cd is a symmetric tensor field that is divergence free and trace 
free with respect to X a b, r is a scalar function, ip is a positive definite scalar function and W c is a 
vector field. Substituting the field decompositions ([9])- (fT0]) into the vacuum constraint equations 
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([I])-© and using standard conformal transformation formulas for the scalar curvature and for 
divergences, we obtain 

(11) (L X W) C 

(12) A A V> 

(Here A\ denotes the Laplace-Beltrami operator with respect to the metric A.) 

The idea of the conformal method is to choose any conformal data (M, A, v, r) in which u is 
traceless and divergence- free with respect to A, and then use the coupled PDE system (fTT|) ~ (fT2]) to 
solve for the determined data (?/>, W). If, for a given set of conformal data, one can solve the system 
(|lip -(|12 p . then the initial data fields obtained by recomposing the fields as in ([9j)- (|10p provide a 
solution to the constraint equations. 

To execute such a construction of initial data one needs to have a symmetric, traceless, and 
divergence-free tensor v. There is a standard method for finding such a tensor [23]. The idea 
behind the method is to start with an arbitrary traceless symmetric 2-tensor field p and then find 
a vector field X which satisfies 

(13) L X X = (div A /i)». 

Using © one easily verifies that v := p, + *D\X is symmetric, traceless and divergence-free. 

There have been extensive studies to determine which sets of conformal data lead to solutions, 
and which do not. (See [4] for a recent review.) This issue is best understood for conformal 
data with constant r, which leads to initial data with constant mean curvature ("CMC"). In 
this case, the constraint equations (fTT]) - (fT2|) effectively decouple — the (unique) solution to ([IT]) is 
T>\W = — and one need only analyze the solvability of the (remaining) Lichnerowicz equation 

(14) A A V - gi?(Aty + \\vcd\lr 7 ~ = 0. 

It is also fairly well understood for "near CMC" conformal data sets, which are characterized by 
|Vr| A being sufficiently small. 

2.3. Conformal method and Andersson-Chrusciel initial data. The Andersson-Chrusciel 
construction [1] of polyhomogeneous AH CMC initial data on a smoothly conformally compact 
3-dimensional manifold M starts by choosing a smoothly conformally compact metric A on M and 
a traceless, symmetric 2-tensor p = p~ l p for which p G C°°(M). The construction continues 
by finding a solution X of the vector Laplacian equation (|13p and by considering the symmetric, 
traceless, divergence-free tensor v := p + T>\X. The polyhomogeneity of X, and consequently of v, 
arises naturally here as a consequence of the indicial roots of the vector Laplacian (for details on 
indicial roots see [18]). More precisely, we have the asymptotic expansions 

(15) X ~ p 2 X + (p 4 logp)Xi with X ,Xt G C°°(M), 

(16) v ~ p^Vo + (plogppi with V ,V 1 G C°°(M). 

Note that the description in [1] treats v as a contravariant 2-tensor, which accounts for the difference 
between our powers of p and the powers of 0, and x in [1]. It is useful that Andersson-Chrusciel [1] 
also prove a sequence of existence and uniqueness results regarding solutions of equations such as 
(1131) in the context of polyhomogeneous tensor fields; we rely on these results when we conclude 
that the perturbations we make are polyhomogenous. 



= \r{\)^ - hu cd + 2{V x W) cd \\^- 7 + i-rV 5 - 
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The construction in [1] proceeds with the analysis of the Lichnerowicz equation. It is shown 
that the boundary value problem consisting of the Lichnerowicz equation (|14p with r = 3 and the 
boundary condition = 1 has a polyhomogeneous solution of the form 

oo 

(17) i> ~ l+p% + Y,(P 3 lQ g PYA, with^GC°°(M). 

i=l 

The exponent in p 3 log p corresponds to the indicial root of the linearization of the Lichenrowicz 
operator (cf. the left-hand side of (|14p with r = 3) in the neighborhood of the constant function 
ipo = 1. We point out that the analysis in [1] also includes an existence and uniqueness result for 
the Lichnerowicz boundary value problem with polyhomogeneous data. We need this result when 
we show that our solution of the Lichnerowicz equation is polyhomogeneous (see Theorem I25|) . 

For convenience, we say that a polyhomogeneous AH initial data set (M, g, K) is of Andersson- 
Chrusciel (A~C) type if it can be written in the form g = and K = g + ip~ 2 v, in which A is a 
smoothly conformally compact metric, v is a symmetric 2-tensor field that is divergence free and 
trace free with respect to A and has an asymptotic expansion of the form (|16p . and ip is a positive 
function with an asymptotic expansion of the form (|17p . The discussion in [1, Appendix A] shows 
that generically, initial data of A-C type are the "smoothest possible" AH initial data. 

Throughout this paper, we assume only that the asymptotically hyperbolic initial data sets we 
work with are poly homogenous in the sense of Definition Q] (which includes data of A-C type as a 
special case). Our gluing procedure then produces new data of the same type. We note that if our 
starting data set is of A-C type, it will not generally follow from our main theorem (Theorem [2] 
below) that the solution we obtain after gluing is also of A-C type; our results will only guarantee 
that this solution is polyhomogeneous in the sense of Definition [TJ The difficulty is that while A-C 
data is obtained by solving the conformal constraints with smooth conformal data, in carrying out 
the gluing we must solve these equations for conformal data which includes log terms. 

3. Main Gluing Theorem 

With the conventions established above, we are ready to state our main theorem. For convenience 
of exposition, we focus our attention here on the case dimM = 3. We have little doubt that the 
theorem and its proof generalize easily to higher dimensions. 

Theorem 2. Let (M,g,K) be a polyhomogeneous asymptotically hyperbolic CMC initial data set 
(with dimM = 3) that satisfies the Einstein (vacuum) constraint equations, and letp%,p2 be distinct 
points in the ideal boundary dM. Then for each e > there exists a polyhomogeneous AH initial 
data set (M e , g e , K £ ) such that 

(i) M e is diffeomorphic to the interior of a boundary connected sum, obtained from M by 
excising small half-balls B\ around p\ and B2 around P2, and identifying their boundaries. 

(ii) (M e , g £ , K £ ) is a solution to the vacuum constraints. 

(iii) On the complement of any fixed small half-balls surrounding p\ andp2 in M , and away from 
the corresponding neck region in M £ , the data (g £ ,K £ ) converge uniformly in C 2 ' a x C ,a 
to (g, K), for some a £ (0, 1). 

In fact, the convergence of K £ is a little better than C > a : away from the fixed half-balls, it 
actually converges in a weighted C 1,a space. See Theorem 1191 for the precise statement. 

Note that, as is the case for the non-localized gluing of AH data sets at interior points (see [15]), 
there is no need to impose any nondegeneracy conditions on the data in the neighborhood of the 
gluing points p\ and p2- 
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4. Splicing Construction 

We presume that we are given a 3-dimensional CMC polyhomogeneous asymptotically hyperbolic 
initial data set (M, g, K), which is not assumed to be connected. We let p denote a chosen smooth 
defining function for M. We may write 

K = g + \i = p~ 2 g + p~ x % 

where g and JL are polyhomogeneous, g G C 2 (M), pZ G C (M), and p, is trace-free with respect to g 
(or, equivalently, g). Note that the assumption that ~g is polyhomogeneous and of class C 2 means 
that the first log term in the expansion of ~g must occur with a power of p strictly greater than 2, 
and thus g is actually in C 2 ' a (M) for some a G (0, 1); and similarly ~fi G C l,a {M). 

The gluing construction is a step by step procedure. We outline the main steps here: The first 
step, which we call splicing, involves the construction of a one-parameter family of manifolds and 
initial data sets that are CMC and polyhomogeneous AH, but that only approximately solve the 
constraint equations. (In most of the literature discussing gluing constructions, both the first step 
leading to approximate solutions and the complete construction leading to exact solutions are called 
"gluing." Here, to distinguish the two, we call the procedure leading to the approximate solutions 
"splicing.") 

The new manifolds M £ are obtained by a connected sum construction which is executed in the 
preferred background coordinates. The parameter e labels the coordinate "size" of the "bridge" , or 
gluing region. Next we use cutoff functions tied to the parametrized gluing region to construct a 
parametrized set of metrics g E on M £ . To verify that these spliced metrics are all asymptotically 
hyperbolic, we also construct a parametrized set of spliced defining functions p £ . Using a different 
cutoff procedure, we produce a family of (spliced) symmetric 2-tensors p £ that, by construction, are 
trace free with respect to the corresponding g £ , but are generally not divergence free with respect 
to g £ . The next two steps involve deformations of the spliced data sets (M £ , g £ , K £ = g £ + p £ ) to 
produce the glued data sets which satisfy the constraints and have the desired limit properties. To 
deform p £ , we first estimate its divergence, and then (following the standard York prescription) 
solve a linear elliptic system (based on the vector Laplacian L 9e ) whose solution tensor deforms p £ 
to a new family of tensors u £ that are divergence free. To deform the metric, we treat (M £ ,g £ , v £ , 3) 
as a set of CMC conformal data, and proceed to solve the Lichnerowicz equation (|14p for a family 
of conformal factors ip £ . The e-parametrized data sets [M £ ,i\> £ g £ ^~ 2 v £ + tp^g £ ), which we call the 
glued data, then solve the constraint equations, and are verified to approach arbitrarily close (as 
e —* 0) to the original data away from the gluing region. 

In the rest of this section, we detail the splicing constructions. We detail the deformation steps 
in subsequent sections. 

4.1. Preferred background coordinates. We focus first on the given polyhomogeneous asymp- 
totically hyperbolic geometry (M,g). Let p be a smooth defining function, and define g = p 2 g, 
which is a C 2 ' a polyhomogenous Riemannian metric on M. For each point p G <9M, we can choose 
smooth functions O 1 ^ 2 such that (p, O 1 ^ 2 ) form smooth coordinates in a neighborhood U C M, 
which we call background coordinates. Sometimes for reasons of notational symmetry we also set 
9° = p. Throughout this paper, we will index such background coordinates with indices named 
a, b, c, . . . , which we understand to run from to 2; and we will use indices j,k, . . . , running from 
1 to 2, to refer to coordinates on dM. We will use the Einstein summation convention when 
convenient. 

It is shown in [12] that when g is an asymptotically hyperbolic metric that is smoothly confor- 
mally compact, there is a smooth defining function p such that \dp\ 2 /p 2 = 1 in a neighborhood of 
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the ideal boundary dM, so the metric can be written in the form g = p 2 {dp 2 + h(p)\ there, where 
h(p) is a smoothly-varying family of metrics on dM with /i(0) = ~9~\ T Qjj- Unfortunately, that result 
does not apply in the present circumstances because we are not assuming that g has a smooth 
conformal compactification. As a substitute, however, we have the following lemma: 

Lemma 3. // (M , g) is a polyhomogeneous asymptotically hyperbolic Riemannian geometry, then 
there exists a smooth defining function p such that 



(18) ^L = 1 + (p 2 ). 

p 2 

Also, for eachp £ dM there exist smooth background coordinates {p,9 l ,9 2 ) on an open neighborhood 
U of p in M in which g can be written in the form 

(19) g = p- 2 [dp 2 + (dO 1 ) 2 + (d9 2 ) 2 + m ab (p, 9)d9 a d9 b 

where the "error terms" m a b are uniformly bounded in IA and satisfy 

(20) m 00 (p,9) = m j0 {p,9) = m 0j {p,9) = 0(p 2 ), je{l,2}, 

(21) m jk (p,O) = O(p+(0 1 ) 2 + (9 2 ) 2 ), j, fee {1,2}. 

Moreover, g is uniformly equivalent inlA to the metric p~ 2 (dp 2 + (d9 1 ) 2 + (d# 2 ) 2 ). 

Proof. Let po be any smooth defining function for M, and write g = p\g. The hypothesis implies 
that there is a smooth metric g^ on M such that g = g 1 + 0(p 2 ). (Just take g 1 locally to be 
equal to the leading smooth terms in an asymptotic expansion for g , and then patch together 
with a partition of unity.) Let gi = Pg 2 5i- Because g\ is asymptotically hyperbolic and smoothly 
conformally compact, the argument of [12] shows that there is a smooth defining function p such 
that \dp\ 2 g x j p 2 = 1. It follows that \dp\g/p 2 = 1 + 0(p 2 ). Let g = i*g = t*~g~i be the metric induced 
on dM by inclusion. 

Given p G dM, let (9 1 , 9 2 ) be Riemannian normal coordinates for g on some neighborhood of p 
in dM. Extend (9 1 ,9 2 ) to a neighborhood of p in M by declaring them to be constant along the 
integral curves of the smooth vector field grad^ p. It follows that (p, 9 1 ,9 2 ) are smooth coordinates 
in a neighborhood U of p, in which has an expression of the form 

g x = dp 2 + (d9 1 ) 2 + (d9 2 ) 2 + m ab (p, 9)d9 a d9 b , 

with moo, m j0, and moj identically zero, and m^ satisfying (f2T|) . Because g = p~ 2 {g 1 + 0(p 2 )), 
this implies that g has the expansion claimed in the statement of the lemma. Since m a fc(0, 0) = 0, 
by shrinking IA we may also ensure that the coefficients m a b are uniformly small in U, and thus g 
is uniformly equivalent to p~ 2 {dp 2 + (d9 1 ) 2 + (d9 2 ) 2 ) there. □ 

From now on, we assume p is a smooth defining function satisfying (|18|) . We now argue that we 
can choose p so that it also satisfies 

(22) A g p < on M. 

On any asymptotically hyperbolic 3-manifold, an easy computation (see, for example, [12, p. 199]) 
shows that 

^£--1 as p^O. 
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Thus, there is some 5 > such that A g p < on the set where p < S. Let a: [0, +oo) — > [0, jS] be 
any smooth, increasing, concave-down function for which 

3 

a(x) = x if x < 5/2, a(x) = —5 if x > 5. 

Define p := a o p. Note that the conclusions of the previous lemma still hold if the function p is 
replaced by p. Furthermore, we compute 

A g p = div s ((cr' o p)dp) = [a' o p) A g p + (a" o p) \dp\ 2 g < 0, 

where the last inequality follows from the facts that a' > 0, a" < 0, and A g p < on the support 
of a' o p. From now on, we replace p by p, and assume that (|22[) holds. 

We call any coordinates {p,9 l ,9 2 ) that satisfy the conclusions of the previous lemma and (|22p 
preferred background coordinates centered at p. 

4.2. Splicing the manifolds and the metrics. We now focus on the topological aspect of our 
gluing construction. Let p\,p 2 & dM be two distinct points on the ideal boundary, and for i = 1,2 
let 0i = (p,9j,9f) = (9®,9},9 2 ) be preferred background coordinates on a neighborhood U C M 
centered at pi. There is a positive constant c such that these preferred coordinates are defined and 
(|19l )-( l2TT) hold for \0i\ < c; after multiplying p and each of the coordinate functions 9\ by 1/c (which 
does not affect (fLSl) . (fT9l) or (f22j) ). we may assume that these two preferred coordinate charts are 
defined for \9i\ < 1. 

We now let e be a small positive parameter, and consider two "semi-annular" regions A e> \, A E £ C 
M characterized by 

A £>i := G Ik : e 2 < |^| < l} . 

We let A e j = ^4 £i j n M. For each choice of e, the two regions can be identified using an inversion 
map with respect to a circle of radius e, given explicitly in coordinates by 92 = I £ {9\), where 
I £ : A e< i — > A e fi is the following diffeomorphism: 

(23) Wi) = -§-B x . 

Fir 

Based on this map, we define an equivalence relation on M by saying 9\ ~ ^2 when 9\ E ^4 e ,i, 
^2 S ^4e,2) and ^2 = ^e(^i)- This produces the connected sum manifold M £ , defined as follows: 

Definition 4. For i = 1,2 and a > 0, Zet B a ^ be the closed subset of M that corresponds in 
coordinates to the ball \9i\ < a. We define Q e C M to be the open subset 

U £ :=I\(B eV UB^), 

and define the spliced manifold M e by 

M £ := 

We let £l £ = £l £ C)M, and let M £ denote the subset of M £ consisting of points whose representatives 
are in £l £ . Let ir £ : £l £ — » M £ be the natural quotient map, and define the neck of M £ to be the open 
subset 

H £ := 7T £ (l e> l) = tt £ (A £t2 ) . 

We let M £ denote Jf £ <T\M £ . 
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We will parametrize the neck by an expanding family of half-annuli in the upper half-space. Let 
H 3 denote the closed upper half-space, defined by 

M 3 :={(y,x\x 2 )eR 3 :y>0}, 

and let M 3 C H 3 be the subset where y > 0. Let r = (y 2 + (x 1 ) 2 + (x 2 ) 2 ) 1 / 2 , and define A £ C H 3 
to be the half-annulus defined by 

(y,^ 1 ,^ 2 ) G H 3 : £ < r < - 



A, 



and let A £ = i £ fl i 3 . Analogously to the case of background coordinates, on H 3 we use the 
notations x° = y and x = (y, x) = (y, x 1 , x 2 ) = (x°, x 1 , x 2 ). 

To define the parametrization of the neck, we first define diffeomorphisms a £: i : A £ — ► A £j i for 
i = 1, 2 by 0i = a £ti (y, x), where 

a £ ,i{y,x) = (ey,ex). 

Then we define (3 £ : A £ —>■ A £t 2 by 

(5 £ = I £ o a £t i = q £i2 ° I, 

where /: A £ — > »4 e is the inversion in the unit circle: I(y,x) = (y/r 2 ,x/r 2 ). Our preferred 
parametrization of jV e is 

^ e := 7T e o a £)1 = TT £ o /3 £ = 7T e o Q ej2 o / : A e -> A4- 

The various diffeomorphisms are summarized in the following commutative diagram: 



«e,l 




The topology of M £ does not change with (sufficiently small) e. The Riemannian geometry on 
M £ (which we define next) does depend on e; this is one of the reasons that we keep track of the 
parameter e. 

To obtain a suitable Riemannian metric g £ on M £ , we blend the metrics coming from the original 
annuli with the use of a cutoff function. 

Lemma 5. There exists a nonnegative and monotonically increasing smooth cutoff function ip: R — ► 
R t/iai is identically 1 on [2, oo), is supported in (^,oo), and satisfies the condition 



(24) 



<^(r) + if 



1. 



Proof. Let <£>o be a nonnegative and decreasing smooth cutoff function such that (fo(r) = | for 
r < \ and 9?o( r ) = for r > 2, and set 



"fir) g ~ ^o(^) + 9^0 
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An easy computation shows that ip satisfies the conclusions of the lemma. □ 

Using this cutoff function and the maps a Sj i, (3 £ , and ^ e defined above, we define the metric on 
M £ as follows. 

Definition 6. We define g £ to be the metric on M £ that agrees with {ir £ )*g away from the neckM £ , 
while on J\f £ it satisfies 

(25) K9e = viDictetfg + tpfy (f3 £ )*g. 

Note the following: 

• On the set where r > 2, ^*g £ agrees with ct* e l g; 

• On the set where r < |, ^>*g £ agrees with (3*g. 

It is obvious from the definition that g £ is poly homogeneous and conformally compact of class C 2 . 

4.3. Splicing the defining functions. Next we construct a family of defining functions for the 
manifolds M £ that are specially adapted to the metrics g £ , and that agree with the original defining 
function p away from J\f £ . To define them, we need the following auxiliary function. 

Lemma 7. There exist a constant A > 2 and a C°° function F : (0, oo) — > (0, oo) that satisfies 

(26) F^=r 2 F(r), r€(0,oo), 

(27) F(r) = 1, r > A, 

(28) F (r) = ^, r<l/A, 
, , \3rF'(r) +r 2 F"(r)\ 1 

(29) l - 4^ ^<^, r G (0,+oo). 

Proof. We use the function 1 + 1/r 2 to interpolate between 1/r 2 for small r and the constant 
function 1 for large r. More specifically, let ipo be any smooth cut-off function such that 

^o(r) = 1 for r < 1, < ^o(r) < 1 for 1 < r < 2, ip (r) = for r > 2, 

let 

C= sup |3r^o(r) + r 2 <(r)|, 

re(0,oo) 

and choose A > 2 large enough that 

The function ip{r) := vpo(Ar) satisfies ifj(r) = 1 for r < 1/A, ip(r) = for r > 2/A and 

(31) |3ry/ (r) + rV'MI < C for all r. 

Define 



F(r):=l + l-^r)-lv 



and compute 



/• ( ij = 1 + r 2 -V ) -r 2 V>(r) = r 2 F(, ■). 
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One readily verifies ([27]). ([28)) and 

( 1 + r~ 2 - if>(r) if re [1/A,2/A] 



(32) F(r) 



Note that 



l + r~ 2 if r e [2/A, A/2] 

l + r" 2 -^V(i) if re [A/2, A]. 



(33) F(r) > if r < 1, F(r) > 1 if r > 1. 



Furthermore, we see that 3rF'(r) + r 2 F"{r) is nonzero only on the intervals [1/A, 2/A] and [A/2, A]. 
A straightforward computation using (|32[) and (|33|) shows 

V 2 |3n//(r) + rV(0l if r € [1/A, 2/A] 

|3r.r (r) + ri 1 (r) 



F(r) 



< < 



13^(1/0 + ^(1/01 if re [A/2, A]. 
Property (|29[) is now immediate from (|30[) and (|3ip . □ 

Note that the function F of the preceding lemma is bounded below by 1, and satisfies an estimate 
of the form 

, \ 1-^(01 c 

v ' F(r) ~ r 

We now let £ : H 3 ->Rbe the positive function defined by 

€(y,x) = yF{r), 

where F is the function of the preceding lemma and r = \x\. A computation using (|26|) shows that 
£ o I = £, and e£ agrees with ey = (a £j i)*p where r > 2 and with ey o I = (fl £ )*p where r < 1/2. 
Define /9 e : M £ — > R to be the defining function that is equal to p away from the neck, and on the 
neck satisfies 

% Pe = e£ = eyF. 

We wish to show that g £ is asymptotically hyperbolic. To do so, we first need to find an expression 
for g £ along the lines of (Tl9|) . Recall that g := y~ 2 {dy 2 + (dx 1 ) 2 + (dx 2 ) 2 ) is the metric of the upper 
half-space model of hyperbolic space. The inversion I(y,x) = (y/r 2 ,x/r 2 ) is an isometry for this 
model: I*g = g. For i = 1, 2, we can write 

dx a dx b 

(35) (a £ji )*g = g + m ab< i(ey,ex) , 

y y 

where m^i are the error terms from the expression (|19p for g in preferred background coordinates 
centered at pj. This metric is uniformly equivalent on *4 £ to g; and because m,^ 6 C 2,a {U n M) 
and 771^(0,0) = 0, it is immediate that as e — > 0, the C 2,a norm of the functions ra ab ^{ey^ex) is 
O(e) on any subset of »4 e where r is bounded above, including in particular the set where ip ^ 1. 
(The C 2 ' Q norm in use here is the ordinary Euclidean one inherited from M 3 .) 
To analyze {(3 £ )* g = I*{a £ ^)* g, we compute, for a = 0, 1, 2, 
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and therefore 

(36) (p £ )*g = I*(a Et2 y g = g+ V m a6>2 (% ^) Q ac (y, x)— — . 

a,b,c,d 

Because Q ac is a rational function with nonvanishing denominator, and is homogeneous of degree 
zero, it is uniformly bounded on A £ , and thus ((3 £ )*g is uniformly equivalent on A £ to g. Moreover, 
all of the derivatives of Q ac are uniformly bounded on any subset of A £ where r is bounded below 
by a positive constant, such as the set where (p ^ 0. Also, on any such subset, an easy argument 
shows that the C 2 '°-norm of rug}, 2{ey / r 2 , ex / r 2 ) is 0(e) as e — > 0. Combining these observations 
with formula (|35p for {a £) 2)*g, we conclude that the pullback of g e to A e has the form 

(37) (V £ )*g £ = g + k ab , £ (y,x) , 

y y 

where for each e, the function k a {, £ is bounded and in C 2 ' a (A £ ), (^ £ )*g £ is uniformly equivalent to 
g independently of e, and 

(38) \\k abi£ \\ C 2, a{Ac) = O(e) 

as e — ► for any fixed c G (0, 1). 

To obtain estimates for the behavior of k ab)£ at the ideal boundary analogous to (|20p and (|2ip . 
we need to explicitly expand the various terms in (|55|) . Note that in addition to the functions Q ab 
being uniformly bounded, we have 

Qi (y,x) = QV(y,x) = ^ = o(^ 

Therefore, using (|20]) - (|2T|) for m abt i, we have 

ho,s{y,x) = cp(r)m o,i{£y,ex) + ip(-) (m jkt2 (^f,^) Q° J {y,x)Q 0k (y,x) 

£y £x\ r .Qj, i \,-,oo/„, „\ i ™ ex \ r>00'„, ^\/o00 



(39) + 2m i0j2 (if, ~ j Q^(y,x)Q 00 (y,x) + m 00 ,2 (^f , =f) Q 00 (y,x)Q 00 (y,x 

,2 s 



^(r)0( £ V) + ^(-)0 



where the implied constants on the right are uniform in e on all of A £ . Note that, in view of the 
definition of p £ , the right-hand side of this equation is 0(p £ ), uniformly in e. Applying similar 
computations to the other terms, and using the notation g abt£ to denote the components of ^* £ g £ in 
standard coordinates on A £ , we conclude that 

(40) gaoAv>*) =y-\i + 0{p £ ))- 

(41) g0j;£ {y,x)=y- 2 O(p £ ), jG{l,2}; 

(42) g jk , £ (y, x) = y- 2 (5 jk + 0(e)), j, k G {1, 2}. 
It follows that the inverse matrix satisfies 

(43) g E °(y,x) = y 2 (l + O(p £ ))- 

(44) g° £ i(y,x) = y 2 0(p £ ), j € {1,2}; 

(45) gi k (y, x) = y 2 {V k + O(e)), j, k G {1, 2}. 
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Lemma 8. There exists a constant C > independent of e such that sufficiently close to the ideal 
boundary dM £ of M £ we have 



\dp. 



Pi 



1 



< —Pe- 

£ 



Proof. Away from the neck, p £ and g £ match the original p and g on M, and the result follows there 
from Lemma [3l We compute on the neck by identifying it with A £ by means of the diffeomorphism 



^ £ . We obtain 



(46) 



\dPe\j 
Pi 



\eF{r)dy + eyF'(r) dr\ g ^ 
(eyF(r)y 

dy 2 , J** nr) dr 

v „, \ y ' F ( r ) 



+ 



F'(r) 



F(r) 



dr 



9e V » v I I 9e 

Using (|43p . we see that the first term on the right-hand side of (|46[) is 

2 

= l + 0(p £ ). 



dy 
11 



To estimate the other terms, note that p e and g £ agree on the neck with (pullbacks of) p and g 
except on the subset Ai/ 2 C A £ . So in the computation below, we may assume that 1/2 < r < 2, 
which means that p £ is bounded above and below by constant multiples of ey, and ^>* £ g £ is uniformly 
equivalent to g. Thus, up to a constant multiple, the second term in (f4"6|) is bounded on Aim by 



dy F'(r) 



F(r) 



dr 



\F'(r)\ C 

77? / \ o — Pe: 

b (r) r z e 



and the third by 



The result follows. 



F'(r) 



F(r) 



dr 



F'(r) 2 2 C 
b (rr e 



□ 



The previous lemma implies that \dp e \ p 2 ge = 1 on the ideal boundary dM £ . Thus, the manifold 
(M e ,g e ) is asymptotically hyperbolic. A consequence of this property is that for each value of e, 
R(g £ ) approaches —6 as p £ — > 0. We need to show that the convergence is uniform in e. 

It follows from [12] that any AH metric g satisfies 

(47) R(g) = -6\dp\ 2 s + p ■ Pi(g,g-\dg) + p 2 ■ P 2 (g,g-\dg, d 2 g), 

where (in accord with the definition of asymptotic hyperbolicity) g = p 2 g and where P m , m = 1,2, 
are certain universal polynomials whose terms involve m -order derivatives of the components of 
g. Since (by Lemma[3|) we have that \dp\^ = 1 + 0(p 2 ), it follows that R(g) + 6 = O(p). Clearly 
then, for each individual £ the function p~ l (R(g £ ) + 6) is bounded on M £ . The uniformity question 
now is whether the C°' a -norms of the functions are bounded uniformly in e. 



c°' a (Ac) ™ bounded uniformly in e. 



Lemma 9. If c > is fixed, then 
Proof. First observe that 

(R(g £ ) + 6)) o xfr E = (eyF)- 1 (R(V* e g e ) + 6) . 



(Ps 1 
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To estimate the scalar curvature R^^gg) let h £ := y 2 ^*g £ . We then have 

(48) R [%g £ ) + 6 = (6 - 6 \dy\\^j + y ■ Px(h £ , hj\dh e ) + y 2 ■ P 2 (h £ , h~\dh £ ,d 2 h £ ). 

The components of the metric h £ can be expressed as 

h a b,e(y, x 1 , x 2 ) = 5 ab + k ab)S (y, x) 

where ||&a&,e||c* 2 > a (A:) = 0(e) as e — > 0. Consequently, the components of h £ , h~ , e~ l dh £ and 
e~ 1 d 2 h £ have uniformly bounded C 0,ct (^4 c )-norms. It follows that the terms in (eyF)^ 1 (R(ty*g e )+6) 
coming from the last two terms of the right-hand side of (|48p have uniformly bounded C 0,a (A c )- 
norms. A careful consideration of the expansion (|39|) yields 



on A c . Therefore, 



koo, £ {y,x) = 0{ey 2 ) and dk 00>e (y,x) = 0(ey) 



e-V 1 (6-6|dy|| e ) 



' X y 1 k 00 , £ (y,x)\\ c0>a{Ac) 



C°> a (A c ) " 

is also uniformly bounded. This observation completes our proof. 



□ 



We conclude the discussion of the AH geometries (M £ ,g £ ) and their defining functions p £ by 
proving that each p £ is superharmonic. 

Lemma 10. If e > is small enough then A gs p £ < 0. 

Proof. Away from the gluing region ty £ {Ai/\), the quotient map tt £ is a diffeomorphism satisfying 
TT £ g £ = g and ir*p £ = p. Thus, away from the gluing region the inequality we need to show is an 
immediate consequence of ([22jh 

To prove the inequality on the gluing region we utilize the transformation law for the conformal 
Laplacian (see, e.g., [19, eq. (2.7)]) to write the Laplace operator for {^ £ )*g £ in terms of that of 
the conformally related metric h £ := y 2 (^/ £ )*g £ : 



3/2 



, v A^N = A,,,, (yF) = l V_ 
Pe J yF 8 F 



^ h Ay l/2 F)- l -R(h £ )y^ 2 F 



It follows from (|38p that the difference between h £ and the Euclidean metric 5 approaches zero, in 
the sense that \\h £ — <^||c 2 > a (.4i/ A ) ~^ 0- Thus we have R(h £ ) — ► and 



,,3/2 



A, (yV*F)-A s (y^ F ) 



V2, 



as e — > 0, 



with both convergences uniform on A.\/\. A straightforward computation shows that 



A 5 (y 1 / 2 F) = --y-^F + 3 

4 r 



-y- 3/2 F- 



t!l F > + v V2 F ». 



Since Lemma [9] shows that R ((^ £ )*g £ ) = —6 + 0(^*p £ ) , which is equal to —6 + 0(e) on Ai/a, it 
follows that 

' A 9ePe 



\]jr* 



Pe 



y 2 3rF' + r 2 F" 
+ 1 — — • 







uniformly on A\/\ as e — > 0. Our result is now an immediate consequence of (|29p and the fact that 
y 2 < r 2 everywhere. □ 
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4.4. Splicing the traceless part of the second fundamental form. Recall that our given 
second fundamental form on M can be written K = fi + g, where p is a traceless, divergence-free 
symmetric 2-tensor field of the form p = p~~ l ~p for some ~p £ C 1,a (M). Our goal in this section is to 
create on each M £ a traceless symmetric 2-tensor p e that is "approximately divergence- free," and 
such that 7r*/x e is equal to p away from the neck. Later we will correct it so that it is divergence- free. 

Let x : R — ► R be a smooth nonnegative function such that x( r ) = 1 f° r r > 3 and supp(x) C 
(2, oo). For each e > 0, define Xe ■ M — > R by 



(49) x £ (P,^) = xf4 + E ML h on ^ 



1,1) 



1,2, 



Then let ju e = XeP on (The level sets of Xe are half-ellipsoids with radius proportional to e 
in the p direction and to yfe in the ^-directions. We have designed this unusual cutoff function so 
that the divergence of ju e will be uniformly small, despite the fact that the tangential and normal 
components of p vanish at different rates near the ideal boundary; see Lemma [TBI for details.) 

It follows from our choice of Xe that ju £ is supported in the set M \ (-B2e,i U -E>2e,2) • Because 
7r e restricts to a diffeomorphism from this set to an open subset of M £ , we can define a symmetric 
2-tensor p £ on M £ by 

(50) p £ = n £ *Ps, 

understood to be zero on the neck. Because g e = ir £ *g on the support of p £ , and \x e is a scalar 
multiple of 7r £!t/ u, it follows that p £ is traceless with respect to g £ . Although it is generaly not 
divergence- free, we will show below that its divergence is not too large (see Lemma [T8|) . 

5. Analysis on the Spliced Manifolds 

In this section we develop the results we will need about linear elliptic operators on our spliced 
manifolds. 

5.1. Weighted Holder spaces and linear differential operators. To carry out the needed 
analysis on AH geometries with AH data, and also to provide a convenient framework for specifying 
the rate at which various quantities like the trace free part of K approach their requisite asymptotic 
values, it is convenient to work with weighted Holder spaces. Here we recall the definition of these 
spaces, using the conventions of [18]. 

Suppose (M, g) is an asymptotically hyperbolic Riemannian geometry of class C^ ,/3 and p is a 
smooth defining function. (Our poly homogeneous metrics, for example, are automatically asymp- 
totically hyperbolic of class C 2 ' a for every a £ (0, 1).) Let (p, , 9 2 ) be background coordinates on 
an open subset IA C M, which we may assume extend to a neighborhood of the closure of IA in M. 
Let B C H 3 be a fixed precompact ball containing (1, 0, 0). A Mobius chart for M (or more accu- 
rately a Mobius parametrization) is a diffeomorphism <3? : B — > U whose coordinate representation 
has the form 

(p, e 1 , 6 2 ) = x) = (ay, ax 1 + b 1 , ax 2 + b 2 ) 

for some constants (a,^ 1 ,^ 2 ). There is a neighborhood W of dM in M covered by finitely many 
background charts, and then the resulting family of Mobius charts covers WDM. We extend this 
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cover to all of M by choosing finitely many interior charts, which we also call Mobius charts for 
uniformity, to cover M \ W. 

Let £ be a tensor bundle over M. For any nonnegative integer k and real number a E [0, 1] 
such that k + a < £ + (3, we define the intrinsic Holder space C k,a (M; E) as the set of sections u of 
E whose coefficients are locally of class C k,a , and for which the following norm is finite: 

|M|fc,« = SU P ||$ ^H^fc.a^); 

where the supremum is over our collection of Mobius charts, and the norm on the right-hand side 
is the usual Euclidean Holder norm of the components of a tensor field on B. For any real number 
8, we define the corresponding weighted Holder space by 



y 5 

with norm 



Cs' a (M; E) = {p 5 u : u £ C k > a (M; E)}, 



\ u \\k,a,6 = \\P S u\\k,a- 



When the tensor bundle is clear from the context, we will usually abbreviate the notation by writing 
C?' a (M) instead of C?' a (M; E). The index 5 labels the rate of asymptotic decay of a given quantity, 
measured in terms of the intrinsic (asymptotically hyperbolic) Riemannian metric g. In particular, 
we note that larger positive values of 5 imply more rapid decay. It is shown in [18, Lemma 3.7] 
that if 77 is any covariant r-tensor field on M with coefficients in background coordinates that are 
(jk,a U p £ Q jjjgaj boundary, then n S Cr' a (M). Similarly, any vector field with coefficients that 
are C k,a up to the ideal boundary lies in Cj"(M). 

A linear partial differential operator P of order m between tensor bundles is said to be geometric 
if the components of Pu in any coordinates can be expressed as linear functions of the components 
of u and their covariant derivatives of order at most m, with coefficients that are constant-coefficient 
polynomials in the dimension, the components of g, their partial derivatives, and 1 / det g% j , such 
that the coefficient of each jth derivative of u involves at most the first m — j derivatives of g. The 
operators A g (the Laplace-Beltrami operator), div 9 (the divergence), T> g (the conformal Killing 
operator), V* (the adjoint of 2? 9 ), and L g (the vector Laplacian) introduced above are all examples 
of geometric operators. It is shown by Mazzeo [20] (see also [18]) that every geometric operator P 
of order m on an asymptotically hyperbolic Riemannian geometry of class C 1 ^ defines a bounded 
linear map from C*' Q (M) to <7™ Q (M): 

(51) \\Pu\\k-m,a,6 < C\\u\\ k>a)S , 

whenever m<k + a<l-\-f3] and moreover, if P is also elliptic, then it satisfies the following 
elliptic estimate for < a < 1 and m<k + a<£ + (3: 

(52) IMka,<5 < C(\\Pu\\ k -m, a ,S + |M|o,O,<0- 

Because our spliced manifolds (M e ,g e ) are poly homogeneous and asymptotically hyperbolic of 
class C 2 , they are also asymptotically hyperbolic of class C 2,a for some a £ (0,1), and thus 
the results we have just discussed hold on M e for each e, with (£,/3) = (2, a). However, for 
our subsequent analysis, we need to check that the constants in (|51|) and (|52p can be chosen 
independently of e when e is sufficiently small. Threading through the arguments of [18], we see 
that for a given geometric operator P, the constants depend only on uniform bounds of the following 
type as $ ranges over a collection of Mobius charts covering M: 

(53) \\<S>*g-g\\ 02 , a{B) <C, sup (^g)' 1 g < C. 

B 
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(In [18], attention is restricted to a countable, uniformly locally finite family of Mobius charts, but 
that additional restriction is used only for Sobolev estimates, which do not concern us here.) In fact, 
it is not necessary to use Mobius charts per se, in which the first background coordinate is exactly 
equal to p\ the arguments of [18] show that it is sufficient to use any family of parametrizations 
satisfying (|53p . as long as there is a precompact subset -Bo of B such that the images of the 
restrictions still cover M, and the following uniform estimates hold in addition to ([53]): 

(54) \\$*P\\c*,<*(B) ^ CPo, ^Po < < C Po , 

where po = < &* / o(l, 0, 0). Thus to obtain our uniform estimates, we need only exhibit a family of 
charts for each e such that the corresponding estimates hold for g = g £ and p = p E , with constants 
independent of e. 

Start with the family of all Mobius charts for M. On the portion of M £ away from the neck, 
these same charts (composed with ir £ ; see Definition [3]) serve as charts for M £ , which satisfy (|53|) 
and (|54p uniformly in e. Recall that we use the diffeomorphism \I/ £ : A £ — > N £ to parametrize the 
neck. Because ^*g £ is isometric to g except on a subset of Ax/i C A £ , we need only show how to 
construct appropriate charts covering points in ^ e (^4 1 / 2 )- 

On this set, we will use standard coordinates on Ai/% as a substitute for background coordinates. 

Given p = ty e (ya, Xq, Xq) £ ^eOAi/2), we define <& p = ty e o tp p : B — ► M £ , where ip p : B — > A £ is the 
map 

<Pp(x, y) = (yoy, yox 1 + xq, y x 2 + xq). 

Note that the Jacobian of ip p is yo times the identity. Under this map, ([37|) shows that g £ pulls 
back to 

K9e =9 + k a b,e(yoy> VOX 1 + Xq, y X 2 + Xq)— — , 

y y 

where we recall that g is the metric of the upper half space model of hyperbolic space. If we assume 
that e is small enough that <p p (B) C A1/4, these metrics satisfy the estimates in (I53p uniformly in 
e because the functions k a t, i£ are uniformly small in C 2,a norm on A\u. The defining function p £ 
pulls back to 

%p £ {x) = £y yF(\<p p (x)\), 

and 3>*p e (l,0, 0) = eyoF(\ (yo, Xq, Xq)\) . Because F is uniformly bounded above and below on A1/4 
by positive constants, and all of its derivatives are uniformly bounded there, it follows that the 
functions & p p £ satisfy the estimates in (|54p uniformly in e. 
Summarizing the discussion above, we have proved 

Proposition 11. Suppose P is a geometric operator of order m < 2 acting on sections of a tensor 
bundle E — > M , and for each e > 0, P £ is the corresponding operator on M £ . There exists a constant 
C independent of e such that for all C 2,a sections u of E, all integers k such that m < k < 2, and 
all real numbers 5, 

H-Pe^llfe— m,a,S — C||^||fc,o,<5' 

If in addition P is elliptic, then 

\\u\\k,a,8 < C(\\P £ u\\k- m ,a,5 + |M|o,0,<s)- 
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5.2. The Vector Laplacian on Hyperbolic Space. In this section, we study the kernel of the 
vector Laplacian L§ = Tft o on hyperbolic space (M 3 ,g). We denote the standard coordinates 

by (y,x) = (y^x 1 ^ 2 ) = (x ,^ 1 ,^ 2 ) on H 3 , and we use the notations \x\ = ((x 1 ) 2 + (x 2 ) 2 ) 1 / 2 and 
r = \x\ = (|x| 2 + y 2 ) 1 / 2 . As a global defining function on H 3 , we use 

ft!/.* 1 .* 2 ) = | x |2 + (y + 1 )2- 

The function p is the pullback to the upper half-space of the usual defining function |(1 — \x\ 2 ) on 
the unit ball. 

It is well known (see, for example, [15] or [18]) that the vector Laplacian 

L s : Cf a {M 3 ) -► C°' Q (H 3 ) 
is invertible for — 1 < 5 < 3 and < a < 1. This leads to the following lemma. 

Lemma 12. If — 1 < 5 < 3 ; i/ien i/iere is no nonzero global vector field X on M 3 satisfying both 
L^X = and the estimate \X\^ < C/5 5 . 

Proof. The hypothesis implies that X 6 C^' (HI 3 ), and then Lemma 4.8(b) of [18] implies that 
X € Cfi ,a (W 3 ) for < a < 1. The result then follows from the injectivity of L§ on the latter 
space. □ 



We need some variations on this result, in which the defining function p is replaced by other 
weight functions. As in Section \A. 31 let £: H 3 — ► R be the function £(y,x) = yF(r), where F is the 
function of Lemma [3 We noted earlier that £ o I = £, where /: H 3 — ► H 3 is the ^-isometry given 
by inversion with respect to the unit hemisphere. Away from and do, £ is a defining function for 
cffl 3 , but it blows up at both and oo. 

Proposition 13. If — 1 < 8 < 3, then there is no nonzero global vector field X on H 3 satisfying 
both L)jX = and the estimate \X\§ < C£ 5 . 

Proof. Suppose X is a nonzero vector field on H 3 satisfying L^X = and \X\§ < C^ 5 for some 
— 1 < 5 < 3. As a consequence of the behavior of the metric g near the ideal boundary, the 
components of X in standard coordinates on H 3 satisfy the condition \X a \ < Cy£ s . 

Let <p: R 2 — > R be a smooth bump function supported in the set where \x\ < \ and satisfying 
< <p(x) < 1, and define a smooth vector field Y = (Y°, Y 1 , Y 2 ) on H 3 by 

(55) Y a (y,x)= [ X a (y,x-u)cp(u)du. 

Differentiation under the integral sign shows that LgY = 0. 

We show first that on the set where r < 1, Y satisfies an estimate of the form \Y\s < Cy s for 
some s with — 1 < s < 3. Observe that by definition £(y,x) = 0(y/r 2 ) for r < 1, and consequently 

\X a (y,x-u)\<Cy( _ V 2 ) for r = \(y, x)\ < 1. 
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Making the substitution u = x — yv (thereby defining v), we have 

8 



ir ' i< - c LA 

-of 

J\x 

= Cy*- S [ 



y 



x — u\ 2 + y 2 

y 



y 



\x-w\<$ \y*\vr+y 

>,- * I dv 

\x-yv\<\ (\v\ 2 + l) 6 ' 



du 



y 2 dv 



Because \(y, x)\ < 1, the triangle inequality implies that {v : \x — yv\ < \} is contained in {v : \v\ < 
2/y}. We now distinguish three cases. 

Case 1: If 5 > 1, then the integrand above has finite integral over all of R 2 . Therefore, 
\Y a \ < Cy 3 - S , from which it follows that \Y\ S < Cy 2 - & . 

Case 2: If 5 < 1, then we let (t,u) denote polar coordinates in the (v l ,v 2 ) plane, and we 
compute 



\Y a \ < y 



3-5 



2tt 



2 /y tdtduj 



o (t 2 + 1)* 

CyZS [(*2 + l)-*+ll^ 2/ " 



<C'y3- 5 (l + y 25 - 2 ) 
< C"y l+S . 



t=o 



It follows that \Y\ S < C"y 5 . 

Case 3: If 8 = 1, then computing in polar coordinates as before, we get 

\Y a \ = Cy 2 [log(t 2 + 1)^ V 
<C'y 2 \logy\<C'y 1+s , 

for any s such that < s < 1. It follows that \Y\§ < Cy s . 

In the three cases above, on the set where r < 1, we have obtained an estimate of the form 
|^|§ < Cy s for some s such that —1 < s < min{<5, 2 — 5}. On the other hand, if r > 1 and \u\ < ^, 
then we have \(y,x — u)\ ~ \(y, x)\ and £(y,x — u) ~ £(y,x), where ~ means "bounded above and 
below by constant multiples of." It follows easily that \Y a \ < Cy^ 5 , and therefore \Y\jj < Ct? on 
this set. 

Now let Y be the vector field Y = I*Y. Because / is an isometry and ^ is /-invariant, the 
argument above implies that 

'Ci\ r<l, 
C(yoI) s , r>l. 



\Y\i< 



Defining a new vector field Z on H 3 by 

Z a (y,x) = [ Y a (y,x-u)<p(u)du, 

JR2 
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we find that L^Z = 0, and consequently the same argument as above shows that Z satisfies the 
estimate 



\z\ s < 



Cy s , r < 1, 

C(yoI) s , r>l 

< C'p s . 

As a consequence of Lemma [12l this implies that Z = 0. 

If X ^ 0, choose a point (yo,xo) S H 3 at which some coordinate component X a (yo,xo) is 
nonzero. After a translation in the x- variables (which is an isometry of H 3 ), we may assume that 
xq = 0. There is some ball B r (0) C M 2 such that X a (yo,x) does not change sign for x £ B r (0). 
If (p is chosen to be supported in this ball, it follows from ([55]) that Y a (yo,0) ^ 0. Repeating this 
argument with Y in place of X shows that there is a point at which Z / 0. This is a contradiction, 
so we conclude that X = as claimed. □ 

We also need the following consequence of this result, in which the weight function is taken to 
be the vertical coordinate y. 



Corollary 14. If — 1 < S < 3, then there is no nonzero global vector field X on H 3 satisfying both 
L^X = and the estimate \X\^ < Cy s . 



Proof. If S > 0, this follows from the previous proposition and the fact that y < C^. If 5 < 0, then 
it follows from Lemma [12] and the fact that y > Cp. □ 

5.3. The Vector Laplacian on the Spliced Manifolds. We now consider the vector Laplacians 
on our spliced manifolds. For each e > 0, let (M e ,g e ) be the asymptotically hyperbolic spliced 
manifold defined in Definitions S] and [6J and let L £ := L 9e be its corresponding vector Laplacian. 
Since g £ is asymptotically hyperbolic of class C 2 ' a for some a £ (0,1), the analysis in [15] or [18] 
shows that L £ : C^' a (M e ) — > C s ,a (M e ) is invertible so long as e > and — 1 < 5 < 3. We need to 
show that the norm of its inverse is bounded uniformly in e. The main goal of this section is to 
understand this uniformity. 

Fix a as above and 5 S (—1,3). We start with the uniform Schauder estimate (see Proposition 

EU) 

( 56 ) ||^||2,a,5 < C'(||^£^||o,a,5 + ||^||o,0,(S) , 

where C is some constant independent of e. We will show that there is a uniform constant D such 
that for sufficiently small e 

(57) ||X|| 0A< 5 < D\\L e X\\ 0fitB . 

This last estimate implies that ||-X"||2 a 5 < C(D + l)||L e X||o a s', i-e., that the norm of the inverse 
(L,)- 1 is bounded above by 1/(C{D + 1)). 

We use blow-up analysis to prove (|57p . The main ingredient in the analysis is the following 
lemma. 

Lemma 15. Let (£,7) be an asymptotically hyperbolic manifold, and let {Nj} be a sequence of 
open subsets of £ such that every compact subset of £ is contained in Nj for all but finitely many 
j. Suppose that for each j we are given a Riemannian metric gj on Nj such that gj — > 7 uniformly 
with two derivatives on every compact subset of £. Assume furthermore that there exist vector 
fields Yj on Nj, a positive real-valued function £ on £, a compact subset Kq C f]j Nj, and positive 
constants C\, C2 such that 



(a) M{C 5 \Yj\ gj ) >C: 

-<5| 
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(b) sup (C a |*iU <c i; 

(c) sup(c-*|L & .^| aj ) -0. 

Then there exist a C°° vector field YonY, and a constant C3 for which 

L^Y = 0, |y| 7 < C3C 6 , Y^O. 

Proof. Let K C X be a precompact open set, and let be a slightly larger precompact open set 
containing K. Since the metrics gj converge uniformly on K (with two derivatives) to 7, we may 
assume that the following estimates hold on K when j is sufficiently large: 

for some constant C independent of j. The function £ is bounded above and below by positive 
constants on K, so it follows that \\Yj \\ H o, P ( K 7 ) and H-^YjH^o.p^ 7 ) ar e bounded uniformly in j. 
Sobolev estimates now imply that 

l|ijf||i?2,p(i<r )7 ) < Ck, 
for some new constant Ck depending on K but independent of j. 

By the Rellich Lemma there exists a subsequence Yj n k of Yj that converges in H 1,P (K,^). For 
p > 3, we have a Sobolev embedding iT ,P (K, 7) — > C '°(-fC, 7). This means that there exists a 
pointwise limit 

Yrt := lim Yj K 

where Y^ £ C°'°(K,7). Note that by construction |Y^| 7 < 2CiC 5 on K. 
Consider a nested sequence of precompact open sets whose union is E: 

K x cKi C K 2 CK 2 C K 3 CK 3 . . . . 

We may use the process outlined above to inductively construct sequences Yj nt K m for each K m , 
such that the sequence Yj nj K m is a subsequence of Yj nt K m _ 1 that converges uniformly on K m . The 
diagonal sequence Yj K n converges uniformly on every compact subset of £ to a continuous limit 
Y on £ that satisfies |Y| 7 < 2C±^ '. The assumption (jaj) ensures that Y ^ 0. So, it remains to 
show that LjY = 0. 

Let X be a compactly supported test vector field on S. Since gj n converges to 7 uniformly on 
suppX with two derivatives, we have 

f (L 7 X,Y) dV 7 = .lim / (L gjn X,Y jn>Kn ) dVy 

Thus y is a weak solution to L 7 Y = 0, and it follows from elliptic regularity that Y £ C°°(S) and 
L 7 Y = 0. □ 

We now focus on verifying inequality (|57p . 



Lemma 16. If —1 < 5 < 3, then there exists a constant D such that for sufficiently small e and 

2 

for all vector fields X £ C s ' (M £ ) we have 

(58) ||-X"||o,o,<5 — II -^e^ II 0,0,5- 
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Proof. Suppose not: Then there exist positive numbers £j — ► and vector fields Xj G C^'°(M £j ) 
such that 

1 1 Xj 1 1 0,0,5 = 1 and \\L Ej Xj || ,o,5 — * 0. 

The fact that ||Xj||o,o,5 = 1 means in particular that for each j there exists a point q £ . G M £j . such 
that 

\Xj{q ej )\ ge .>\p Sj {q ej ) 5 . 

Since 7r e maps Jl e / 3 = M \ (-Be/3,1 U -B e / 3i2 ) surjectively onto M e , for each j we may choose a 
representative gj for g £j such that qj G r2 e / 3 . Passing to a subsequence if necessary, we may assume 
that qj — > g G M. Our proof now splits into several cases depending on the location of q. Each 
case culminates in a contradiction. 

Case 1 : 5 £ M. This is the easiest of the cases as it allows immediate use of Lemma [T5l 
Indeed, let (£,7) = (M,g), Nj = fi e ./3> = C = P> an d let i^o be a compact set containing 
a small neighborhood of q. Vector fields Yj on iVj for which {is £ .\Yj = Xj necessarily satisfy the 
hypotheses of Lemma [15j Thus there exists a C°° nonzero vector field Y on M with L g Y = and 
\Y\ g = 0(p s ). However, for 5 G (—1, 3), the vector Laplacian has no kernel in C?'°(M), so this is a 
contradiction. 

Case 2: q G <9M \ {pi,P2}- Let := (p, 9 1 ,9 2 ) be a set of preferred background coordinates 
centered at (7, which we may assume to be defined on a half-disk D C M whose coordinate radius 
is R. Let (pj,6j,0j) be the coordinates of qj, j 3> 0. Note that (pj,#j,#|) — ► (0,0,0) as j — > 00. 

Let (£,7) be the hyperbolic space (M. 3 ,g). We define Nj to be the half-ball in H 3 centered at 
(0, —9j j pj, —9j j pj) of (Euclidean) radius R/pj- As soon as j is large enough that \9j\ < R/2, the 
triangle inequality shows that the half-ball of radius R/2pj centered at (0, 0, 0) is contained in Nj, 
so we see that (J^ Nj = £ and that each compact subset of M is contained in Nj for all but finitely 
many j. 

Consider the transformations Tj : Nj — » M whose coordinate representations are given by 
(59) Tj(y, x\x 2 ) = ( Pj y, Pj x l + 9}, P jX 2 + 9 2 ). 

These transformations are chosen so that 0, 0) = qj. We will now construct metrics and vector 
fields on Nj satisfying the hypothesis of Lemma [151 

First, let gj := T?g. It follows easily from Lemma [3] that gj — > g uniformly on compact sets 
together with two derivatives. 

Now define Yj := p~ s T*Xj. We compute: 

\YjW\g, = Pj S \Xj{^))\ 9 < Pj S (pjV) S = V s , 
1^(1,0,0)1^. = P J S \Xj(qj)\ g > \pfp] = \, 

y~ S Wj\ gj = { P] y)- 5 (\L g X 3 \ g o Tj) < \\L £j Xj\\ 0AS -+ 0. 

In particular, the hypotheses of Lemma [151 are fulfilled for Kq = {(1,0,0)} and £ = y. It follows 
that there is a nonzero vector field Y on H 3 for which LgY = and \Y\g- = 0{y s ). For 5 G (—1, 3) 
this is impossible by Corollary 1141 

Case 3: q G {^1,^2}; without loss of generality we may assume that q = p\. For sufficiently 
large j, the point q £ . is contained in the neck M £j ; let Xj := (yj,Xj,x 2 ) be the point in A £j such 
that ^> £j (xj) = q £j , and let Tj := \xj\. It follows from the fact that qj G £\/3 that Tj > 3. 

There are several different ways in which qj can converge to p\. We consider now three subcases, 
and use Lemma [15] in each subcase. 
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Case 3a: There are uniform upper and lower bounds on r 3 and yjjrj] i.e., for some d > 0, 

1 v- 1 

(60) - <n<d and 1 > ^ > - > 0. 

3 rj a 

In this case we take (£,7) = (H 3 ,g), Nj := .4^ and gj := *I'*.g £j . It is immediate that |J ■ Nj = S, 
and that every compact set in £ is contained in almost all Nj. It follows from (|37p that (fy — > g 
uniformly on compact sets together with two derivatives. 

Consider the function £(y,x) = yF(r) and vector fields Yj such that (^f S j)*Yj = £ J S Xj. Because 

^tjPej = £ j£i we have 



\Yj\ g . < £ 5 and sup (C S Wj\ g ,) -» 0, 

so conditions (jb]) and (jcj) of Lemma 1151 are satisfied with C — £• The compact set i^o characterized 
by (|60p contains the points (x,-), where we have 

\Yj(xj)\ 9j =^l^(^)l 9ej > = 

This means that the condition (jaj) of Lemma [JS] also holds. Therefore, there exists a nonzero 
C°° vector field on Y on H 3 that satisfies LgY = and \Y\g- = O(^). However, this contradicts 
Proposition [TBI 

Case 3b: There is a uniform positive lower bound on yj/rj, but rj are unbounded. Passing to 
a subsequence, we may assume that 

d 



r 3 ■ —> 00, — > - > 0. 



We again take (£,7) = (H 3 ,c/). Consider the transformation Tj : H 3 — > H 3 given by ^j(ic) = r^x. 
This transformation is chosen so that the points dj := T^ 1 (xj) = (yj/rj,Xj/rj,Xj/rj) lie in a 
compact region Kq of the upper hemisphere {r = l,y > 0}. 
The set iVj C H 3 characterized by 



1 1 

Nj := { x : — < r < 



3 r j £ j r j 



is taken via Tj to the outer portion of the expanding annulus A ey Since in this case rj — ► 00 and 
0, we see that (J ■ iVj = H 3 and that any compact subset of HI 3 is contained in almost all Nj. 



3' 3 



Define gj := T*^*.g £j ] because EjTj — > 0, a simple argument using ([37|) shows that gj converges 
uniformly to g on compact subsets of H 3 together with two derivatives. 
Consider 

( = y= £ jrjF(rjr) T ^^ ^ Yj ^ r ^V% X r 
Since F is bounded above and below by positive constants on [1/3, 00), we have 

V- S \Yj\ gi = (ejrjyr S (\Xj\ ge , o %oTj) < ( Ej r jV ) S {T*%p ej ) S = F{r 3 rf < C; 

y~ 5 \L 9j Yj\m = ^rjyr S {\L as .Xj\ ge . o vl; o Tj) < F(rjr) s \\L ej Xj\\ 0AS - 0. 

Moreover, since o Tj{dj) = q £j and p £j (q ej ) = ^Pejixj) = EjyjF(rj), we have 

\Yj(aj)\ 9j = {sjrjT S \Xj{ qEj )\ 9ej > \{e 3 r 3 Y & p £] {q £] ) s = \{y 3 /r 3 ) 5 F(r 3 ) s > Cy{d 3 ) 5 . 
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Consequently, the conditions of Lemma [T5l are fulfilled. This means that we now have a C°° nonzero 
vector field Y on M 3 for which L S Y = and |Y| 3 = 0{y 5 ). This is a contradiction to Corollary [TD 
Case 3c: There is no positive lower bound on yj/rj. Passing to a subsequence, we may assume 
that one of the following holds: 



(1) rj — * oo, or 

(ii) yj — > and rj is bounded. 



Note that in both cases (JTJ) and (jnj), \{Xj,x 2 )\/yj — ► oo. In either case, we consider transformations 
Tj : H 3 -> M 3 defined by 

Tj(y,x x ,x 2 ) = (yjy^yjx 1 + x^^jx 2 + x 2 ), 
chosen so that * E , o7}(l, 0, 0) = g e .. Let iV,- := where .4+ := {r > l}nAe v The region 



3' 



Nj C H 3 is the semiannular region centered at (0, —x^/y 
radius 1/ejyj. 

In case ([1]), we have |(xj,Xj)| — > oo and 



;/%') °f inner radius 1/Qyj and outer 



-xj 



Vj Vj 



1 



2y, 



OO. 



In case ([n|), once j is big enough, we have |(xj,x|)| — | > as a consequence of > g. It follows 
that 



Vj Vj 



1 



(x),x 2 ) 



> 



Vj 



1 

7 %' 



oo. 



In particular, Nj contains the half-ball of radius 



Rq 



1 



mm ■ 



Vj Vj 



centered at the origin. Since EjXj — > 0, the first expression in the minimum also converges to oo. 
Therefore, Rj — ► +oo, |J . iVj = H 3 , and every compact subset of H 3 is contained in almost all Nj. 

As in the previous case, we take (S, 7) = (H 3 ,</) and 5^ := T*^*.g £j . Note that ([57]) shows that 
gj can be expressed as 

9j=9 + k a b,e 3 (Tj(y, x)) — — , 



where fc a b )£ are functions defined on A £ such that \\k, 



ab,e\\c 2 ^(A c ) 



for any fixed c > 0. We need 



to show that \\k, 



ab,s ° Tj\\C 2 < a (K) 



for any fixed compact set K C H 3 . We will consider cases ^ 
and (jn]) separately. Let K C H 3 be a compact set, and let R be the supremum of t\k- 

First assume we are in case flij). For any point x = (y, x) G if, as soon as j is large enough that 
.R < ^| (xj, Xj)\/yj, the reverse triangle inequality gives 

\(x),x 2 )\ 



r oTj (x) 







+ x 



> 



(4j) 



R > 



\yy 



00. 



Eventually, therefore, the set Tj{K) lies in the portion of A e where r > 2, and thus for x G if we 
have 

k a b, £j {Tj{xj) = m abA (ejyjy, EjyjX 1 + ~ '.jUj-''' + £jX 2 ), 
Tj\\ C 2, a (K) ~^ because (£jyj,ejXj,EjX 2 ) — > 0. 



and it follows easily that ||fc, 
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On the other hand, if we are in case ([n]), then for any x G K, as soon as R < \(xj,xj)\/yj, we 
have 

/ Kx 1 x 2 )\ \ 

r o Tj{x) < y 3 f 3 j J ' 1 + R j < 2\(x},x 2 )\ < 2 rj < C, 

since {rj} is bounded. Therefore, Tj(K) is contained in the fixed annulus {x : | < r < C}, on which 
k a b,s converges to zero in C 2,a norm. Because yj — > 0, the transformations Tj are affine transfor- 
mations with uniformly bounded Jacobians, and so again we conclude that ||fc a b,e ° ^j\\c 2 ' a (K) ~* 0- 
In both cases, therefore, gj —* g in C 2 ' a (K). 
This time, we let 

Reasoning as in Case 3b, since i 7 is bounded above on [g,oo) and bounded below everywhere, we 
find that 



y~ 8 



y~ 5 \Yj\ 9j = ! -,//,.</) 'vv ; , o o 25) < ( £i y iy )- 5 (r;**.p £ .) d = f{t o t,) 5 < c ; 

= l v.'/,.'/; o * o T ; ) < Fir o T.Ol^X.-Ho^ - 0; 

1^(1,0, 0)| 9 . = (ejy^lX^)^ > \{e ]yj )- 5 p £3 {q ej f = \yj S { yj F( yj )) S > c. 

These estimates show that the vector fields Yj satisfy the conditions of Lemma [15] for the choice of 
Kq = {(1, 0, 0)}. We now see that there exists a nonzero C°° vector field Y on H 3 such that 

L S Y = 0, Y = 0(y 5 ) for 5 G (-1, 3). 

However, this is impossible by Corollary [Til □ 

Now we are ready for our main theorem concerning the vector Laplacian. 

Theorem 17. If e is sufficiently small and if —1 < 5 < 3, 0<a<l, then the vector Laplacian 

L £ : Cf a (M £ ) -+ C°i a (M £ ) 

is invertible and the norm of its inverse is bounded uniformly in e. 

Proof. The invertibility follows from the analysis in [15] or [18], and the uniform estimate follows 
by combining (|56p with (|58p . □ 



6. Correcting the traceless part of the second fundamental form 

In this section, we use the elliptic PDE theory and analysis of the previous section to add a 
correction to our spliced tensor \x £ to make it divergence-free. First we show that its divergence is 
not too large. 

Lemma 18. With [i £ defined by ([50]) . div Se p £ G C®' a (M £ ) with norm \\ div 9e /U £ ||o,ai,l = 0{y/e) . 

Proof. Recall that we have defined ju £ = XeM, where \i is the given traceless second fundamental 
form and \e is defined by ([49]) . Restricted to the support of ju e , the projection ir £ is a diffeomorphism 
taking g to g £ and % to fx £ , so it suffices to show that div 9 ju e G C 1 ' a (M) with 0{y/e) norm. 

For a vector field Y and a symmetric 2-tensor 77, let us use the notation Y j 7/ to denote the 
1-form rj(Y,-). It is easy to check (by doing the computation in Mobius coordinates) that the 
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map (Y,rj) i— ► Y jr/ is a continuous bilinear map from C^ a (M) x Cg^ a (M) to C^'" 52 (M) for any 
5 h 5 2 G M. 

It follows easily from the definition of the divergence operator and the fact that p is divergence- 
free that 

div 9 % = Xe div 9 p + (grad 9 x e ) J /x = (grad 9 x e ) J M- 
The support of grad g (Xe) is contained in the union of the two half-balls B\ t i U B\p- Letting 9 3 
denote either 9\ or 9\ depending on which half-ball we are in, we compute 

and therefore 

(61) div g % = X ' f^ + E^-j ^(grad gP )j / x + ^^(grad g ^)j^ 

Using the formula for the change in Christoffel symbols under a conformal change in metric (see, 
for example, [18, equation (3.10)]), we find that 

= div 9 p = p 2 divjp - p(gr&dgp) j p. 

After substituting p = p" 1 ^, this becomes 

= p 2 divg(p _1 p) - p(gradgp) J (p _1 p) = pdiv^ p - 2(grad^p) jp. 

It follows that (grad^-p) Jp = ^pdivg p, and thus 

(grad^ p) j p = {p 2 gradg p) J (p~ l ~p) 

= \p 2 dlv-gp. 

Since div^p is a 1-form whose coefficients in background coordinates are in C 0,a (M), div^-p is 
contained in C® ,a (M), and thus (grad ff p) J p G (^'"(M). On the other hand, a straightforward 

computation shows that grad fl 9 J G ^'"(M), and therefore (grad^ 9 3 ) J p G (^'"(M). We conclude 
that the following quantities are finite: 

||p~ 2 (grad 9 p) J/i|| 0l a,i, ||p _1 (grad s ^) J//|| 0lQi) i (j = 1,2). 

With this in mind, we rewrite (1611) as 



(62) divgp E = / e (p,0) j ^-p 2 (grad ff p) jp + ^^pp x (grad ff J ) Jp 



where 



fe{P,0)=X + Z^^ 



Note that / e is bounded independently of e, and is supported in a region where p < e\/3 and 
1^1 < V^£- Its differential satisfies 

„ J p 2 x ^ (6^) 2 \ / 2pdp x - 20W 
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Because \dp\ g and \d9^\ g are both bounded by multiples of p, it follows that \df £ \ g is bounded 
uniformly in e. Therefore, f £ is uniformly bounded in C 1 (M) and thus also in C ' a (M). Inserting 
these estimates into ()62p . we find that 

II div g p s \\ 0l a,i 

< 11/^110^,0^ — ^ ||^ Vad g p)j^|| cw + ^ - \\p H&^Mo,^) 

< C'{e + v/i) 

This completes the proof. □ 

Using the preceding result and Theorem 1 171 the idea now is to make a small perturbation of p e , 
which we denote by u e , for which div Se z^ e = 0. 

Theorem 19. For each sufficiently small e > 0, there is a poly homogenous symmetric 2-tensor 
field v £ , which is traceless and divergence-free with respect to g £ , such that p 2 v £ has a C 2 extension 
to M £ that vanishes on dM £ , and v £ satisfies 

(63) \v E - Aielli.ci = 0{y/e). 

In particular, away from the neck, v £ converges uniformly in C\ ,a {and therefore also in C 1,a ) to 
(the projection of) p. 

Proof. We use the standard technique for finding a divergence-free perturbation of p £ discussed in 
Section [2] (see (fl3j) ). Relying on Theorem [T71 for each small e > 0, we let X £ be the unique vector 
field in C 2 ' a (M £ ) that satisfies 

(64) L £ X £ = (div gE p £ )K 

and we set v £ := p £ + T> £ X £ . By definition of the conformal Killing operator, v £ is traceless; and 
by construction it is divergence- free. Since ||div ffe /x e ||o,a,i = 0(y/e), it follows from the uniform 
estimate of Theorem 1171 that ||X £ ||2 )a ,i = 0(y/e). The arguments of Section 15.11 show that T> £ is 
bounded from C 2,a (M £ ) to C^ ,a (M £ ) uniformly in e, and therefore (|63|) is satisfied. It now remains 
to show that v £ is poly homogeneous and that p £ u £ has a C 2 extension to M £ which vanishes on the 
ideal boundary. 

We start by observing that the right hand side of (|64p is polyhomogeneous and that, by Theorem 
6.3.10 of [1], there exists a polyhomogeneous solution Xf hg of (|64p . Since p has an asymptotic 
expansion beginning with and the first log term (if any) appearing with p s , s > 0, a computation 
shows that the vector field on the right-hand side of (|64p has an expansion beginning with a p 2 
term, and with the first log terms (if any) appearing in the p £ +3 term, s > 0. Inserting the general 
asymptotic expansion for Xf h& into (|64p and matching like terms inductively, we conclude that 
Xf hg has an asymptotic expansion beginning with p 2 and the first log terms appearing with p £ +3 , 
s > 0. (Note that the first log terms which arise from the indicial roots of L £ appear with pf.) On 
the other hand, X £ £ C 1 ' a (M £ ) implies that its component functions in background coordinates are 
also 0(p 2 ). The uniqueness part of Theorem 6.3.10 in [1] implies that X £ = X<P hg . It now follows 
easily that v £ is polyhomogeneous and that is has an asymptotic expansion starting with p~ l and 
the first log term appearing with p|, s > 0. Thus the extension of p 2 £ v £ to M £ is actually of class 
C 2 and vanishes on the ideal boundary, which is just what is needed for v £ to be the traceless part 
of the second fundamental form for a polyhomogeneous AH initial data set (see Definition [1]) . □ 
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7. The Lichnerowicz Equation and Conformal Deformation to the Spliced 
Solutions of the Constraint Equations 

Thus far, starting with a set of asymptotically hyperbolic, polyhomogeoneous, constant mean 
curvature initial data (M, g, K) satisfying the Einstein constraint equations, together with a pair of 
points {P11P2} both contained in the ideal boundary dM, we have first produced a one-parameter 
family of spliced data sets (M £ , g £ , fj, £ ) which are asymptotically hyperbolic, poly homogeneous, 
CMC, and not solutions of the constraints, and we have then corrected fj, £ to a new family of 
symmetric tensors v e which are all divergence free as well as trace free with respect to g £ . We 
have verified that outside of the gluing region, v £ approaches the original trace free part of K in 
an appropriate sense. 

To complete our gluing construction, we will now carry out a one-parameter family of conformal 
deformations that transform the data (M £ ,g £ , u £ , r = 3) to a family of data sets (M e , ^ £ g £ ,ip~ 2 v £ + 
iptde) satisfying the desired properties of the gluing construction (including the constraint equa- 
tions) for all e. Following the principles of the conformal method outlined in Section \2.2\ if we want 
the conformally transformed data sets to satisfy the constraints, then the conformal functions ij) e 
must solve the Lichnerowicz equation (I14j) . which for the data (M £ ,g £ ,u £ ,T = 3) takes the form 
C e (t/j e ) = 0, where 

(65) C £ (u) := A 9e u - ^R(g £ )u + \\v £ \ 2 ge vr 7 - \u\ 

Hence, we need to do the following: prove that for each e the Lichnerowicz equation (|65j) does 
admit a positive solution ip £ which is polyhomogeneous and C 2 up to the ideal boundary, prove 
that ip £ approaches 1 at the ideal boundary (so that the resulting Riemannian manifold is AH) and 
prove that as e — > the solutions ip £ approach 1 away from the gluing region. We carry out these 
proofs here. 

The first step in our proof that, for the data sets (M £ , g e ,u e ,T = 3), the Lichnerowicz equation 
admits solutions with the desired asymptotic properties, is to estimate the extent to which the 
constant function tpQ = 1 fails to be a solution of (|65p . While it is relatively straightforward to 
show that C e (ijjo) is an element of the weighted Holder space C 1 ' a (M £ ), we have been unsuccessful 
in proving that the corresponding norm of C e (tpo) is "small". Consequently, we are able to find a 
solution ip £ of the Lichnerowicz equation such that i/j £ — ?/>o "vanishes" on dM £ but we are only able 
to obtain good estimates on ip £ — ip$ in C 2,a (M £ ). This, however, is sufficient to prove our main 
result. 

Lemma 20. We have £ £ (V>o) G C°' a (M £ ) and ||£ e (V>o)||o,a = 0{y/e) as e -> 0. 
Proof. Note that tpQ = \ yields 

Ce{^)=-\{R{ge) + Q-\Ve\l). 

The fact that C £ (t/jo) G C®' a {M £ ) is immediate from \ v £ \ 2 ge £ C°' a (M e ) which is true by construction, 
and (R{g £ ) + 6) £ C^' a (M £ ) which is true by virtue of (|37|) and Lemma [H To estimate the 
unweighted norm of C £ (i(jq), let <I> be one of our preferred charts for M £ (see Section [5.ip . If the 
image of $ is away from the gluing region, i.e., if &(B) n ^ £ (A = 0, then 

||$*£ e (^o)|| c o,a ( B) = gll$*(N;L - \^\l e )\\ c o, a( B) = 0(Ve) 
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as a consequence of the second constraint equation R(g) — \K\g + r 2 = R(g) — \p\g + 6 = and 



(|63p . Thus, it remains to study the charts for which 

where c > is some sufficiently small fixed number. We start with the inequality 

||$*£ e (V?o)|| C o,« ( b) < h\**( R (9e) + 6)|| C o,« ( b) + 0(Ve) + |||$*|m £ |LHco,« (j b)- 

It follows from (I47p and Lemma[9]that \\R(g £ ) + 6||o, a ,i is bounded uniformly in e. The uniformity 
properties (|54h and the fact that p £ < \fsjc on fy £ (^A c ^j imply that 

(66) \\**{R(9e) + 6)[| c o,a ( ^ ) = 0(Vi). 

To understand the /i e -term, note that 7r e is a diffeomorphism on the support of ju £ , where we also 
have 

|2 . I |2 21 — 1 2 

It then follows that sup |$*|/i e |^ e | = O(e). Likewise, 

B 

\d(\n 2 g )\ g = P \d{p 2 x 2 mi)\- a < p{Mdp\mi + 2p 2 \dxm 2 a + ? 2 kn|)i§) = <V) 

as a consequence of the boundedness of \dp\g and the fact that p 2 \dx\g = pO(e)0(-) = O(p) on the 
support of dx- Overall, we see that 

\\$*\Pe\l\\ C °,c.(B) =°( £ ) 

and therefore ||<&*£ e (^o)|lccw,B) = 0(y/e) as e — ► 0. □ 
It should also be pointed out that \u £ \g s = R(g £ ) + 6 + 8£ e (ipo) implies that 

(67) \\We\l L <C 

v ' II 1 c 'He 1 10, a ~~ 

for some C > 0. 



The main ingredient in our study of the solvability and the solutions of the Lichnerowicz equation 
is the uniform invertibility of the linearizations 

V £ :=A gE -^R(g £ ) + 7\u £ \l+30 

of C £ at tpo = 1. In what follows we rely heavily on maximum principle(s). Part of the reason why 
this approach is successful is that the function 

1/-, - -. ,2 



f e :=-(R(g £ ) + 71^ + 30, 
has a positive lower bound. 

Lemma 21. Let C < 3 be a positive constant and let e be sufficiently small. We have 

fe > C, pointwise. 
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Proof. It is enough to show 

(68) sup \f £ - \u £ \ 2 - 3| -»• as e -> 

or, equivalently, that the sup-norms of \R{g £ ) — We\g s + 6| both over the gluing region ^ £ (A 
and over its complement M £ \ ^(.A converge to 0. To prove this convergence on ^ £ (A ^/s), 
note that tt £ maps diffeomorphically onto the support of fi e , and that 

1 12 ^ ( 2 1 — \2\ -1 
Thus, there is a constant c > such that \l^e\^ e < ce on 

* e (^ /3 ). In light of dSD this means 

that 

1 1 2 

SUp l^elg ~ y as ^ — ^ 0. 

The convergence result (j68j) on the gluing region now follows from Lemma [9] or, rather, estimate 
(|66p . The convergence away from the gluing region is an easy consequence of the fact that the 
restriction of f £ to M £ \ ^ £ [A satisfies 

(fe ~ We\ 2 ga ) °vr £ = | + 7|(7r e )*i/ e g + 3o) - |(7T e )*I/ e g = 3 + ± - |(7r e )*i/ e |j) = 3 + 0(^) 

by virtue of ([63]) and the second constraint equation R(g) — \K\ 2 + 9 = R(g) — \p\ 2 + 6 = 0. □ 

Strictly speaking, the operators V £ are not "geometric" due to the presence of the |^ £ |g e term, 
so the analysis of [18] does not apply directly. There are many ways to circumvent this; for 
convenience, we will base our argument on Proposition 3.7 of [12]. First we need the following 
uniform estimate (called the "basic estimate" in [12]). This estimate is analogous to the estimate 
of Lemma [TU] above for the vector Laplacian. The proofs of the two lemmas, however, are quite 
different: Lemma [TU] is proved using blow-up analysis, while the proof of the next lemma is direct 
and constructive. Consequently, the next lemma features a more optimal result on C (as opposed 
to the vector Laplacian case where we are only able to prove the existence of C). 

Lemma 22. Let C > | be fixed, and assume e > is sufficiently small. 

(1) If u is a C 2 function on M £ with both p £ 1 u and p~ V £ u bounded then 

(69) sup l/OfT 1- "! < C sup \p~ x V £ u\. 

(2) If u is a bounded C 2 function on M £ with V £ u bounded, then 

(70) sup\u\ < C sup\V £ u\. 

M £ M e 

Similarly, if U is a precompact subset of M £ , and u is a continuous function on Q that is C 2 in $7 
and vanishes on d£l, then 

(71) sup|p~ 1 u| < C sup I /3 ~ 1 T^e ii I , and sup|u| < C sup\V £ u\. 

n n n n 

Proof. We start by proving (|69p . Note that it suffices to consider functions u for which 

sup l/Oj 1 ^! = sup (p~ u) • 

Me Me 

Given a fixed e > and a C 2 -function u G C^'°(JW e ), Yau's Generalized Maximum Principle [12] 
implies that there is a sequence of points {xk} of M £ such that 
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(i) lim [p E 1 u](x k ) = sup [p £ 1 

(ii) lim Ic^p" 1 ?/)! (x k ) = 



u 



k^oc 



9e 



(iii) lim sup A Se [p e 1 u] (x k ) < 0. 



k— >oo 



Note that the condition (ii) can be re-written as 



(72) 



lim 

k— >oo 



-1, -1 d Pe 



(x k ) = 0. 



A short computation shows that 
(73) 



a f —1 1 -1 A 0/ d Ps -1 , -1 d Pe\ -1 A g E Pe 

A 3, LPe u J = Pe A 9e U - 2 \ — >Pe du - Pe It—) - p e U— . 

Pe Pe Pe 



Recall that, by Lemma [HJ the quantity 
identity ([72]) implies 



is bounded for each fixed e > 0. Therefore, the 



(74) 



lim sup 

k— >oo 



-1 A -1 A 9ePe 

Pe A 9e U ~Pe u ~ 

Pe 



(x k ) = lim sup A ge [p £ 1 u] (x k ) < 0. 

k^oo 



Since the defining functions p £ for small e > are superharmonic (Lemma [TO]) and since f £ > ^ 
by Lemma (2TJ we see that 



Pe lA 9e U ~ Pe lu 



Conditions (i) and (|74|) now imply 



p- l V £ U + Pe" 1 



> Pe X V e U + ^Pe V 



lim sup [pe l V £ u] (x k ) + ^ sup [p £ 1 u] < 

fc^+oo M e 

for small enough e. Consequently, we have 

sup [pe -1 ^ < C liminf [ — p~ x V e v\ (x k ) < C sup \p~ 1 V £ u\, 

Me k^+oo Me 

as claimed. The proofs of the remaining three estimates are similar but considerably easier. Indeed, 
to prove (f70j) we use 

Ag £ u(x k ) = V £ u(x k ) + f e {x k )u(x k ) > V e u{x k ) + jju(x k ) 
in place of (173p . while (jTlh is proved using the ordinary maximum principle. 



□ 



Theorem 23. The operators T E : Cl' a {M e ) C° ,Q: (Me) and V e : C 2 ' a (M £ ) -» C7 ' Q (M £ ) are 
invertible for sufficiently small e > 0. T/ie norm of the inverse of V £ : C 2,a (M e ) — ► C°' Q (Me) is 
bounded uniformly in e. 

Proof. It follows from Proposition 3.7 in [12] (together with Lemma [22]) that V £ is invertible when 
e is small enough, so it remains only to prove uniformity of the norm. We start by establishing a 
uniform elliptic estimate 

( 75 ) \\ u h,oL < C (||7Vu||o,a + IM|o,o) 

in which C is independent of (sufficiently small) e > 0. Let be one of our preferred charts for M e 
(see section I5TTT) . Consider the elliptic operator V$ tS : C 2,a (B) — ► C 0,a (B) defined by 



V %£ := A^* ge -$*(3 + |i/ e g e -£e(^o)); 
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this operator is of interest since 

(V £ U) = V^^U. 

Recall that the metric Q*g £ is uniformly equivalent to the hyperbolic metric g. Furthermore, we see 
from Lemma I2U1 and (|67|) that the C 0,Q (-B)-norms of $*|i/ £ |^ e and <&*j£ £ (tpo) are uniformly bounded. 
Thus, the eigenvalues of the principal symbol of "P<j> i£ are uniformly bounded from below, while the 
C°' a (-B)-norms of the coefficients of are uniformly bounded from above. Let Bq be a fixed 
precompact subset of B such that the restrictions of our preferred charts to Bq still cover M £ . It 
follows from the standard elliptic theory [11] that there is a constant C (independent of e, 3> and 
u) such that 

Wu\\ c2 , a{&o) < C (\\V*, s $*u\\c*,«(B) + ll**«llco.o ( ^)) • 

Taking the supremum with respect to $ now yields (|75p . 

Next, we combine Lemma [221 and the elliptic estimate ([75]). We conclude that there is a constant 
C (independent of e) such that 

IM| 2 ,a < C\\V £ u\\ , a 

for all u e C 2 > a (M £ ). This shows that the norm of V~ x : C°' a (M £ ) C 2 > a (M £ ) is bounded 
independently of e. □ 

We solve the Lichnerowicz equation by interpreting it as a fixed point problem. More precisely, 
consider the quadratic error term 

QM := £ £ (V'o + r/)-£ e (Vo)-^ = ^ e | 2 ((l + r/)- 7 -l + 7r ? ) - jj((l + r,f - 1 - 5r?) 

and the corresponding map 

Q £ : V ^-{V £ )- 1 (C £ (^ Q ) + QM). 
Note that, by Lemma [50] and Theorem \T6\ 

G £ : C 2 > a (M £ ) - C 2 ' a (M £ ) and Q £ : C 2 ' a (M £ ) - C 2,a (M £ ). 

It is easy to see that the solutions ip £ = ipQ + % of the Lichnerowicz equation correspond to the 
fixed points r\ £ of Q £ . In what follows we argue that Q £ is a contraction mapping from a small ball 
in C 2 - a (M £ ) to itself. 

Lemma 24. For sufficiently large C and sufficiently small e, the map Q £ is a contraction of the 
closed ball of radius Cyfe around in C 2 ' a (M £ ). 

Proof. Let 771,772 G C 2 ' a (M £ ) be of norm 0(y/e). Assuming in addition that |t/i| < 1 and 1 772 1 < 1, 
using the bound on \^e\^ expressed in (|67|) . and using the binomial expansion formulae, we find 
that for sufficiently small e > 0, 



\Qe(m) - Qe(rn)\\o, a 



3 

'4 L 



1 + m)~ 7 - (1 + my 7 + 7(7/1 - m 

(1 + 7 ? i) 5 -(l + 7 ?2 ) 5 - 5(7/1 -7/ 2 ) 



0,0 



< 0(Ve)||7/2 - 7/l||o, Q < 0(y/e)\\ 

A similar calculation shows that if 1 1 77 1 1 2,ck = 0(y/e) then ||Q e (7/)||o, a = 0(e). As a consequence 
of Lemma [20l functions 77 with [|7/||2a = 0(y/e) also satisfy 

\\C 6 (lh) + Q e (v)\\0,a = O(y/i). 
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Combining this with Theorem [23] we have that there is a sufficiently large constant C > such 
that for sufficiently small e > 

WGMh* = ||(^) _1 (A(i) + Q £ (r/))|| 2)Q < CVi- 

Thus we have determined that Q £ : B c ^ — > B c ^. 
To see that this map is a contraction, we compute 

\\G £ (m) - Ge(m)h,a = ||(^)- 1 (Q e (r?i) - 2 £ (%))|| 2jQ = 0(^)llm - r»[| 2>a . 

It thus follows that if e > is small enough, the map Q £ : B c ^ — > B c ^ is a contraction. □ 

We are now ready to state and prove our main result regarding solutions of the Lichnerowicz 
equation for the parametrized sets of conformal data (M £ ,g £ ,v £ ,T = 3): 

Theorem 25. // e is sufficiently small, there exists a poly homogeneous function ip £ on M £ which 
has a C 2 extension to M £ that is equal to 1 on dM £ , and satisfies 

(76) A 5e Ve - lR(gs)A + ~ \j>\ = 0. 

The function ip £ is a small perturbation of the constant function ipQ = l in the sense that 

HVfe - V'olkc* = 0(y/e) as e -> 0. 
Proof. By the Banach Fixed Point Theorem, the sequence 

V0,e ■= 0, r] lj£ := Ge(r)Q,e), ■■■ ,Vn,e-=Ge(lln-l,e),--- 

converges in B c ^ C C 2,a (M £ ). Thus, there exists a function r\ £ on M £ such that ||?y e ||2,a < Cy/e 
and such that the function ip £ := V'o+^e solves the Lichnerowicz equation. To address the regularity 
of tp £ , note that r\ n ^ £ S C^' Q (M e ) for all n,e. Consequently, each r\ n ^ £ has a continuous extension 
to M £ that vanishes on dM £ . Because convergence in C 2,a (M £ ) implies uniform convergence, it 
follows that, for each fixed e, the limit r/ e := lim n ^oo r} n . £ also has a continuous extension to M £ 
and vanishes on dM £ . We now conclude that i/j £ = tpQ + rj £ approaches 1 at the ideal boundary. 
Therefore, Corollary 7.4.2 of [1] applies and we see that ip £ is polyhomogeneous. Inserting the 
asymptotic expansion for ip E into (1T6|) and comparing like terms inductively, we find that the first 
log terms in ip e appear with p^. (These terms arise as a consequence of the indicial roots of the 
linearized Lichnerowicz operator.) It follows that tp £ has a C 2 extension to M £ . □ 

With these solutions ip £ to the Lichnerowicz equation in hand, we readily verify that the one- 
parameter family of initial data sets (M e , t/j^g £ ^~ 2 i , £ + tlj^g £ ) satisfies the list of properties outlined 
in [2j Hence we have constructed the desired asymptotic gluing of AH initial data satisfying the 
Einstein constraint equations. 

8. Conclusions 

The gluing construction which we have discussed and verified here allows one to take a pair 
(or more) of CMC initial data sets for isolated systems with unique asymptotic regions-either 
asymptotically null data sets in asymptotically flat spacetimes, or data sets in asymptotically 
deSitter spacetimes-and glue them together in such a way that the spacetime which develops 
from this glued data has a single asymptotic region. In the case that the original data sets are 
asymptotically null, one may wonder how the Bondi mass [22] for the glued data compares with 
the Bondi masses for the original data sets. We will study this issue in future work. 

There are a number of ways in which the results proven here might be extended. It should be 
straightforward to be able to handle solutions of the Einstein-Maxwell or Einstein-fluid constraints, 
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rather than the Einstein vacuum constraint equations. A more challenging generalization we plan 
to consider is to allow for initial data sets which do not have constant mean curvature. We have 
done this in earlier gluing work [10] using localized deformations of the original data sets so that, 
in small neighborhoods of the gluing points, the mildly perturbed original initial data sets do have 
constant mean curvature. The work of Bartnik [3] shows that this sort of deformation can always be 
done. A key first step in generalizing our results here to non CMC initial data sets is to generalize 
Bartnik's local CMC deformation results to neighborhoods of asymptotic points in AH initial data 
sets. This issue is under consideration. 

One further generalization of some interest is to attempt to carry out localized gluing at asymp- 
totic points in AH initial data sets. To do this, it would likely be necessary to determine if the 
work of Chrusciel and Delay [8] generalizes so that it holds in asymptotic neighborhoods in AH 
initial data sets. While this may prove to be difficult, we do believe that we will be able to localize 
the gluing to the extent that in regions bounded away from the ideal boundary, the glued data is 
unchanged from the original data. 

References 

[1] Andersson, L., Chrusciel, P. T., Solutions of the constraint equations in general relativity satisfying "hyperboloidal 

boundary conditions" , Dissertationes Math. (Rozprawy Mat.) 355 (1996). 
[2] , On "hyperboloidal" Cauchy data for vacuum Einstein equations and obstructions to smoothness of scri, 

Comm. Math. Phys. 161 (1994), 533-568. 
[3] Bartnik, R., Regularity of variational maximal surfaces, Acta Math. 161 (1988), 145-181. 

[4] Bartnik, R., Isenberg, J., The constraint equations, The Einstein equations and the large scale behavior of 

gravitational fields, 1-38, Birkhauser, Basel, 2004. 
[5] Beig, R., Chrusciel, P. T., Schoen, R., KIDs are non-generic, Ann. Henri Poincare 6 (2005), 155-194. 
[6] Choquet-Bruhat, Y., Theoreme d'existence pour certains systemes d'equations aux derivees partialles non 

lineaires, Acta Math. 88 (1952), 141-225. 
[7] Choquet-Bruhat, Y., York, J. W., Jr., The Cauchy problem, in "General relativity and gravitation," Vol. 1, 

edited by A. Held, Plenum, New York, 1980, 99-172. 
[8] Chrusciel, P. T., Delay, E., On mapping properties of the general relativistic constraints operator in weighted 

function spaces, with applications, Mem. Soc. Math. Fr. (N.S.) 94 (2003). 
[9] Chrusciel, P. T., Delay, E., Lee, J. M., Skinner, D. N., Boundary regularity of conformally compact Einstein 
metrics, J. Differential Geom., 69 (2005), 111-136. 
[10] Chrusciel, P. T., Isenberg, J., Pollack, D., Initial data engineering, Comm. Math. Phys. 257 (2005), 29-42. 
[11] Gilbarg, D., Trudinger, N.S., Elliptic Partial Differential Equations of the Second Order, Springer- Verlag, 1983. 
[12] Graham, C. R., Lee, J. M., Einstein metrics with prescribed conformal infinity on the ball, Adv. in Math. 87 
(1991), 186-225. 

[13] Isenberg J., Constant mean curvature solutions of the Einstein constraint equations on closed manifolds, Class. 

Quantum Grav. 12 (1995), 2249-2274. 
[14] Isenberg, J., Maxwell, D., Pollack, D., A gluing construction for non-vacuum solutions of the Einstein- constraint 

equations, Adv. Theor. Math. Phys. 9 (2005), 129-172. 
[15] Isenberg J., Mazzeo R., Pollack D., Gluing and wormholes for the Einstein constraint equations, Comm. Math. 

Phys. 231 (2002), 529-568. 
[16] , On the topology of vacuum spacetimes, Ann. Henri Poincare 4 (2003), 369-383. 

[17] Isenberg, J., Park, J., Asymptotically hyperbolic non-constant mean curvature solutions of the Einstein constraint 

equations, Class. Quantum Grav. 14 (1997), A189-A201. 
[18] Lee, J. M., Fredholm operators and Einstein metrics on conformally compact manifolds, Mem. Amer. Math. Soc. 

183 (2006). 

[19] Lee, J. M. and Parker, T. H., The Yamabe Problem, Bull. Amer. Math. Soc. 17 (1987). 

[20] Mazzeo, R., Elliptic theory of differential edge operators I, Comm. Partial Differential Equations 16 (1991), 
1615-1664. 

[21] Mazzeo, R., Pacard, F., Maskit combinations of Poincare- Einstein metrics, Adv. Math. 204 (2006), no. 2, 
379-412. 



ASYMPTOTIC GLUING 



37 



[22] Wald, R. M., General Relativity, University of Chicago Press, 1984. 

[23] York, J. W., Jr., Conformally invariant orthogonal decomposition of symmetric tensors on Riemannian manifolds 
and the initial-value problem of general relativity, J. Mathematical Phys., 14 (1973), 456-464. 



